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Analytic QCD models are those versions of QCD in which the running coupling parameter a{Q^) 
has the same analytic properties as the spacelike physical quantities, i.e., no singularities in the 
complex plane except on the timelike semiaxis. In such models, a{Q^) usually differs from its 
perturbative analog by power terms ~ {IS^ /Q^)^ for large momenta, introducing thus nonperturba- 
five terms ~ (A'^/Q^)*' in spacelike physical quantities whose origin is the UV regime. Consequently, 
it contradicts the ITEP operator product expansion philosophy which states that such terms can 
^ come only from the IR regimes. We investigate whether it is possible to construct analytic QCD 

O models which respect the aforementioned ITEP philosophy and, at the same time, reproduce not 

just the high-energy QCD observables, but also the low-energy ones, among them the well-measured 
semihadronic r decay ratio. 
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P-l I. INTRODUCTION 



Today one of the main goals in strong interaction theory is to technically enlarge the applicability of QCD to 
processes involving lower momentum transfer q^. Thereby several obstacles have to be overcome. One of them is 
that the running QCD coupling a{Q^) = as{Q'^)/iT, when calculated within the perturbative ("pt") renormalization 
group formalism (we call it Opt), in the usual ("perturbative") renormaUzation schemes, yields singularities of a^tiQ^) 
at > 0, usually called Landau singularities. Consequently, spacelike observables expressed in terms of powers of 
T-H apt(Q^) obtain singularities on the spaceHke semiaxis < < (Q^ — —q^, with q denoting the typical momentum 
transfer within a given physical process or quantity) . This is not acceptable due to general principles of local quantum 
\^ field theory [T. Furthermore, studies of ghost-gluon vertex and gluon self-energy using Schwinger-Dyson equations 
[2] and large- volume lattice calculations 3J, result in QCD coupling a{Q^) without Landau singularities at > 
and even with a finite value at Q = 0. Consequently, the behavior of the coupling a{Q^) at low values of should 
be corrected relative to that given by perturbative reasoning. 
^ Several attempts at achieving such corrections have been recorded during the last 14 years starting from (what 

we call) the minimal analytic (MA) QCD of Shirkov and Solovtsov [3]. Here, the trick lay in simply omitting the 
wrong (spacelike) part of the branch cut within the dispersion relation formula for a{Q'^). Consequently, the resulting 

^ analytized coupling A^^^\Q^) = a('^^^((3^) is analytic in the whole Euclidean part of the plane except the 
nonpositive semiaxis: G C\(— cx),0]. Furthermore, for evaluation of physical observables which are represented, in 
ordinary perturbation theory, as a (truncated) series of powers of apt(Q^), one also has to extend the analytization 
procedure to a^^. {n > 2). In MA this was performed in Ref. [5] (see also Ref. 0) and resulted in the replacement 

of 

o-pt nonpower expressions A^'^\Q'^). This specific procedure was dubbed by the authors of [SI |S] analytic 
perturbation theory (APT); whereas we will refer to it generally as minimal analytic (MA) QCD. 

Other analytic models for a{Q^) satisfy certain different or additional constraints at low and/or at high [TUlSj. 
Analytic QCD models have been used also in the physics of mesons [El within the Bethe-Salpeter approach, and 
in calculation of analytic analogs of noninteger powers a'^^ jl8' within the MA model (for reviews of various analytic 
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QCD models, and further references, see Refs. [THHH])- We note that the MA couphngs An {n > 1) defined here 
are the MA couphngs of Refs. [H |6l [20] divided by tt. 

Ah of these versions of analytic QCD have one common feature: their (analytized) coupling a{Q^) differs from the 
perturbative coupling even at higher energies by a power term: 

|<5a(Q')I^KQ')-«pt(g')l^(AVQY (Q'»A2), (1) 

where A; is a positive integer (usually = 1; for the models of Refs. [T^I15| : k — 3). How can these power corrections 
be interpreted? In a given (usual) renormalization scheme, where apt(Q^) has (Landau) singularities on the positive 
axis ^ A^(~ 0.1 GeV^) > 0, analytization of apt(Q^) can be understood to be achieved by a modification of 
the discontinuity ("spectral") function Pi^{(t) = Imapt(Q^ = —a — ie) at energies \cr\ < A^, thereby subtracting 
the Landau singularities from apt{Q^). It is this subtraction, in the given renormalization scheme, which leads to 
the power deviations Eq. ([l]) and, as a consequence, to terms ^ {h? /Q'^)^ in all spacelike physical quantities. But 
such contributions are definitely of nonperturbative origin, since they are proportional to e'xp{— K / a-pt{Q'^)) which is 
nonanalytic at Opt = [cf. Eq. (fTol) in Sec. [TT|. 

Whether such terms, producea in spacehke observables 'D{Q^), can be interpreted as being of ultraviolet (UV) 
origin or not, is not entirely clear. Interpretations of such terms in the literature differ from each other. For example, 
Ref. [22] suggests that the Landau pole is not of (entirely) UV origin because the Landau pole persists in the 
renormalization group resummed expression for apt(Q^) even if one uses, instead of UV logs, the mass-dependent 
polarization expression (with a sufficiently small gluon mass). On the other hand, the authors of Ref. [23] argue 
that the aforementioned terms ~ (A^/Q^)'^ are of UV origin due to the following consideration: If one considers the 
leading- /3o summation of an inclusive spacelike observable V{Q'^) (cf. Appendix O) 

V^^^\Q^)^ r^Fj,{t)a{tQ^e^), (2) 
Jo * 

where Fxi{t) is a characteristic function of the observable and C ~ —5/3, then the quantity tQ^e'^ indicates the 
magnitude of the (squares of) internal loop momenta appearing in the resummation. In the UV regime of these 
momenta, e.g., for i > 1 (see also Ref. [21]), the deviation ([I]) then leads to power terms of apparently UV origin in 
the observable 

^ Fvit) ^ (Ayq')' . (3) 

Considering all these arguments, we come to the conclusion that the aforementioned (A^/Q^)*^ contributions in 
physical quantities are at least partially due to UV effects. The existence of nonperturbative contributions stemming 
from the UV regime is not in accordance with the operator product expansion (OPE) philosophy as advocated by 
the ITEP group [23l[25]. This philosophy rests on the assumption that the OPE, which has originally been derived 
in perturbation theory (PT), is valid in general (i.e., even when including the nonperturbative contributions) and 
consequently allows for a separation of short-range from long-range contributions to (inclusive) QCD observables. 
While the short-range contributions can be calculated perturbatively and lead to expressions for the OPE coefficient 
functions, the long-range contributions show up as matrix elements of local operators and can be parametrized in 
terms of condensates (not accessible by PT). And it is this long-range part which leads to power corrections reflecting 
the contributions of nonperturbative origin to the observable. Therefore, according to the ITEP interpretation, the 
power term corrections stem from the IR region. This ITEP-OPE approach rests on intuitive physical arguments, 
and has led to the success of QCD sum rules. 

In this work we will adopt the aforementioned ITEP philosophy when analytizing perturbative QCD and, conse- 
quently, we will request that the analytic coupling parameter Ai{Q^) = a{Q^) differ from the usual perturbative one 
at high by less than any power of A^/Q^. 

We wish to stress, however, that there is nothing in quantum field theory (QFT) that would impose on us the ITEP 
interpretation of the OPE. In this context, we mention that the essential singularity at a = [such as exp(— ii'/a)] 
has quite a general and mysterious genesis - first mentioned in QFT by Dyson |26j on specific physical grounds, and 
later by many authors on more formal grounds (for an overview, see 2T and references therein). 

An additional feature of most versions of analytized QCD is that they fail to reproduce the correct value for the most 
important (since most reliably measured) QCD observable at low energies, namely Tt, the strangeless semihadronic r 
decay ratio, whose present-day experimental value is (cf. AppendixjB]): rT-(exp.) = 0.203 ±0.004. Most of the analytic 
QCD models are either unable to predict unambiguously value, or they predict significantly smaller values (e.g., in 
MA, Ref. [SJI2H]), unless unusual additional assumptions are made, e.g., in MA that the light quark masses are much 
higher than the values of their current masses [21] . 
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This finding (loss in the size of rr) in MA appears to be connected with the eUmination of the unphysical (Euclidean) 
part of the branch cut contribution of perturbative QCD. Since r-r is the most precisely measured inclusive low 
momentum QCD observable, its reproduction in analytic QCD models is of high importance. The apparent failure 
of the MA model with light quark current masses to reproduce the correct value of r,- had even led to the suggestion 
that the analytic QCD should be abandoned [30j. 

Here, we are investigating whether a modified version of QCD can be defined which simultaneously fulfills the 
following requirements: 

(i) It is compatible with all analyticity requirements of Quantum Field Theory. In particular, it must not lead to 
Landau singularities of a(Q^), and furthermore we expect (see Sec. [n]) that a{Q^) is analytic at = 0, and 
thus IR finite, with a{Q^ = 0) = oq < oo. 

(ii) It is in accordance with the ITEP-OPE philosophy which means that the UV behavior of a{Q^) is such that 
\a{Q'^) - apt(Q^)| < [t^ jQ"^)^ for any integer k at large Q^. 

(iii) The theory reproduces the experimental values for Vr (and other low energetic observables, e.g. the Bjorken 
polarized sum rule at low Q^). 

We will show that such a theory is attainable, but only at a certain (acceptable, we think) price. Some of the main 
results of the present work have been presented, in a summarized form, in Ref. |31| . 

We are approaching our aim in an indirect way, namely by properly modifying the /3 function /3(x) \x = a((5^)] 
of QCD. This approach, which has been used first by R§,czka [32 in a somewhat different context, means that the 
starting point in the construction is the beta function /3(a), rather than the coupling parameter a{Q^) itself or its 
discontinuity function pi(cr) = Ima((5^ = —a — ie). The ITEP-OPE condition can be implemented in such an 
approach in a particularly simple way (see below). Consequently, we are trying to augment /3(a) which, in general, is 
only specified by its perturbation series around the point a = 

/3(a) = -fSoa'^il + cm + c-^a^ + cga^' + 0{a% , (4) 

where /3o and c\ — Pi/ [3q are two universal constants. This should be done in such a way that the augmented beta 
function leads (via the renormalization group equation RGE) to an effective analytic coupling a{Q^) which also enables 
the correct evaluation of low-energy QCD observables in a perturbative way. 

The abovementioned requirements for a{Q'^) imply the following constraints on the modified beta-function /3(a): 

(1) The /3 function must be such that the RGE gives a running coupling a{Q'^) analytic in the entire complex plane 
of Q^, with the possible exception of the nonpositive semiaxis: E C\(— oo,0]. 

(2) For small |a|, /3(a) has Taylor expansion Q in powers of a, i.e., the perturbative QCD (pQCD) behavior of 
/3(a), with universal /3o and ci, at high is attained. 

(3) /3(a) is an analytic (holomorphic) function of a at a = in order to ensure |a(Q^) — apt((9^)| < (A^/Q^)''' for 
any A: > at large (see Sec. 0), thus respecting the ITEP-OPE postulate that powerlike corrections can 
only be IR induced. At high Q^Tthose pQCD values apt{Q^) which reproduce the known high-energy QCD 
phenomenology are attained by a{Q^). 

(4) It turns out to be difficult or impossible to achieve analyticity (holomorphy) of a(Q^) in the Euclidean complex 
plane E C\(— oo,0] unless the point = is also included as a point of analyticity of a{Q^). This then 
implies that a{Q^) — >■ ao when — >■ 0, where qq is finite positive, and that /3(a) has Taylor expansion around 
a = ao with Taylor coefficient at the first term being unity: /3(a) = (a — qq) + 0{{a — ao)^). Then, /3(a) is a 
nonsingular unambiguous function of a in the positive interval a E [0,ao]- Note that analyticity of a((5^) at 

= is in full accordance with the general requirement that hadronic transition amplitudes have only the 
singularities which are enforced by unitarity. 

We proceed in this work in the following way. In Sec. |II] we construct various classes of beta functions which give 
analytic a(Q^) at all E C\(-oo,0) and fulfill the ITEP-OPE condition. We relegate to Appendix |a] details of 
the analytic expressions for the implicit solution of RGE and their implications for the (non) analyticity of a{Q'^). In 
Sec. |III| we point out the persistent problem of such models giving too low values of r,-. In Sec. |IV| we present further 
modification of the aforementioned beta functions, such that, in addition, the correct value of r,- is reproduced. In 
Appendix |B] we present the extraction of the massless and strangeless t-t value from experimental data. We relegate 
to Appendixes [C][d] and |e] the presentation of formalisms for the evaluation, in any analytic QCD (anQCD) model, of 
massless observables, such as the Bjorken polarized sum rule (BjPSR), the Adler function and the related r,-. Appendix 
[C] presents construction of the higher order anQCD couplings; Appendix [P] presents a formalism of resummation of 
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the leading-/3o (LB) contributions in anQCD; Appendix [e] presents a calculation of the beyond-the-leading-/3o (bLB) 
contributions in anQCD. Section [V| contains conclusions and outlines prospects for further use of the obtained anQCD 
models. 



II. BETA FUNCTIONS FOR ANALYTIC QCD 



Our starting point will be the construction of certain classes of beta functions /3(a) for the coupling a{Q^) such 
that ITEP-OPE conditions 



|fl(Q')-apt(Q')|< (^) , (fc = l,2,...), (5) 



are fulfilled and that, at the same time, they lead to an analytic QCD (anQCD), i.e., the resulting a{Q'^) is an analytic 
function for all G C\(— oo,0]. This procedure is in contrast to other anQCD models which are usually constructed 
either via a direct construction of a{Q^), or via specification of the discontinuity function Pi(<t) = lma{Q^ — ~a — it) 
and the subsequent application of the dispersion relation to construct a{Q^) 

In such approaches, it appears to be difficult to fulfill the ITEP-OPE conditions ([5|^ , and difficult or impossible to 
extract the beta function /3(a) as a function of a. 

On the other hand, starting with the construction of a beta function /3(a), which appears in the RGE 

it turns out to be simple to fulfill conditions (Is]) (cf. Ref. [M])- Namely, if one requires that /3(a) be an analytic 
function of a at a = 0, then the corresponding a(Q^) respects the ITEP-OPE conditions ([5|. 

This statement can be demonstrated in the following indirect way: assuming that the conditions ^ do not hold, 
we will show that /3(a) must then be nonanalytic at a = 0. In fact, if the conditions ([s]) do not hold, then a positive 
no exists such that 

a(Q2)«apt(Q2) + ^(AVQ2)"« (8) 

for ^ A^. Asymptotic freedom of QCD implies that at such large the perturbative a^tiQ^) has the expansion 
(if the conventional, MS, scale A = A [33l [34] is used) 



1^2. _ 1 _ ci Inln(QVA^) / ln^(ln(gVA^)) \ 

"P'^^^ /3oln(QVA2) Plh^\Q^/K^)^ \ In^Q^K^) J' 

and consequently the power term can be written as 

(AVQ')"" = exp {-K/aptiQ^)) (/3oapt)-^' (l + ^(aln^ a)) , (10) 
where K = no//3o and K' — noCi//3o. Applying d/dlnQ^ to the relation ([8| and using expression (10), we obtain 

/3(a(Q2)) « /3pt(apt(Q2)) _ no^exp {-K/a^tiQ^)) (/3oapt)-^' (l + 0{aln^ a)) . (11) 



Replacing a{Q^) in the first beta function in Eq. (11) by the right-hand side (rhs) of Eq. ([8|, using Eq. (10), and 
Taylor expanding the /3(a((5^)) function around apt(Q^) (7^ 0), gives 

/3(apt,) + «;exp(-i^/apt)(/3oapt)-^' (l + ©(aln^ a)) x ^^L=„^^ + O (^exp{^2K/apt)a-,^''') 

«/3pt(apt)-«oKexp(-i^/apt)(/3oapt)-^' (l + 0(aln2a)) . (12) 



^ Instanton effects can modify the conditions (jsjl in the sense that these conditions remain vahd only for A; = 1, 2, ... , fcmax where 2fcmax is 
the largest dimension of condensates not affected by the small-size instantons. Scenarios of instanton-antiinstanton < 4/3o 

(= 9 for rif =3), cf. Ref. |23| . In this work we do not consider such possible instanton effects. 
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In this relation, valid for small values of |apt|, the term with derivative d(3{a)/da ~ apt on the left-hand side (Ihs) 
can be neglected in comparison with the corresponding term on the rhs. Therefore, Eq. ( 12 1 obtains the form (with 
notation a,. 



^ a) 



I3{a) sa /3pt(a) - noKexp(-ii'/a)(/3oa) 



-K' 



1 



©(aln^ a)) 



(13) 



We note that /3pt(a), being a polynomial, is analytic at a = 0. The term proportional to cxp{—K/a) is nonanalytic 
at a = 0, because exp{—K/a) has an essential singularity there. This shows that nonfulfillment of the ITEP-OPE 
conditions ^ implies nonanalyticity of /3(a) at a = 0, and the demonstration is concluded. 

This proof shows that nonfulfillment of ITEP-OPE conditions implies nonfulfillment of a = analyticity of /3(a). 
Or equivalently, fulfillment of a = analyticity of /3(a) implies fulfillment of the ITEP-OPE conditions ([S]). This does 
not mean the equivalence of a = analyticity of /3(a) with the ITEP-OPE conditions. But that will suffice for our 
purpose, since in the following we will simply restrict the Ansdtze for the /3 function which are analytic at a = 0, thus 
having the ITEP-OPE conditions secured. 

Integration of RGE ^ must be performed for all complex Q^. To achieve this, we first need an initial condition 
[equivalent to the fixing of scale (^ 0.1 GeV^)]. This is a subtle point within our approach, due to two reasons. 
First, when we choose a specific form of the beta function /3(a), we automatically choose a specific renormalization 
scheme (RSch) as well, as represented by the coefficients cj = [ij/ {j > 2) of the power expansion of /3(a), Eq. (IT 



The running of the corresponding a{Q^) can be in general significantly different from the running a{Q^ ;MS) in MS 
RSch. Secondly, this running is also influenced by the number of active quark flavors and by flavor threshold effects. 
In our analyses of RGE with our specific /3 functions, we will consider the number of active quark fiavors to be = 3, 
i.e., the flavors of the three (almost) massless quarks m, d and s. We do not know how to include in a consistent way 
the massive quark degrees (n/ > 4) in anQGD. On the other hand, the ITEP-OPE conditions ([s]) tell us that the 
considered anQGD theories become practically indistinguishable from pQCD at reasonably high energies ^ A^. 
Therefore, we wish to keep n/ = 3 in the RGE running to as high values of \Q'^\ as possible, and to replace the theory 
at higher |Q^| by pQGD, in the RSch dictated by the specific beta function. Furthermore, in pQGD the threshold for 
Uf = 3 1-^ Uf — 4: can be chosen at ^ [knicY with k k, 1 — 2> [551 - I55] . where rric denotes the mass of the charmed 
quark. We will use fc = 3, i.e., at > (3toc)^ (~ 14.5 GeV^) the anQGD theory will be replaced by pQCD theory. 

In order to find the value of a((3TOc)^) = ajn which will define our initial condition, we start from the experimentally 
best known value of the coupling parameter, namely a(M^, MS). It is deduced, within pQGD, from all relevant 
experiments at high \Q'^\ > 10^ GeV^ and found to be a(M|,MS) w 0.119/7r, Ref. We RGE run this value, 

in MS RSch, down to the scale (3mc)^, and incorporate the quark threshold matching conditions at the three-loop 
level according to Ref. [38 at = ^ -y^g obtain^ a = a{{3mc)^ ,MS,nf = 3) = 0.07245. The value 

flin = a((3TOc)^), at the same renormalization scale (RScl) but in the RSch as defined by our /3(a) function, is then 
obtained from the aformcntioned MS value a = a((3mc)^, MS, = 3) by solving numerically the integrated RGE in 
its subtracted form (Ref. ^40^ , Appendix A there) 



1 
a 



- Ci In 



Cia 
1-l-cia 



dx 



x'^{l+cix)(3{x) 



1 



cia 



= + ci In :; ] + dx 

a \l + ciaj Jq 



l3{x) + Pax'^jl + cix) 
x'^{l + cix)]3{x) 



(14) 



where a = a((3mc)^) = a-m and a = a((3mc)^, MS) = 0.07245, both with Uf — 3; further, /S is the beta function of 
the MS scheme. We note that in Eq. (14 1 our beta functions have expansions around a = [cf. Eq. with the 
RSch coefficients (c2, C3, . . .) which may be considerably different from the MS coefficients (c2,C3, . . .). Therefore, in 
Eq. ( 14) expansions of /3 in powers of x are in general not justified. 

(3mc)^) fixed, RGE (l7| can be solved numerically in the Q^-complex 



Having the initial value ain = (^{Qt 



plane. It turns out that the numerical integration can be performed more efliciently and elegantly if, instead of Q^, a 
new complex variable is introduced: z = ln((5^//ifn)- Then the entire Q^-complex plane (the first sheet) corresponds 



to the semiopen stripe — tt < Im^; < -\-tt in the complex z plane. The Euclidean part Q G 



-00, 0] where a{Q^) 



has to be analytic corresponds to the open stripe — tt < Im(z) < +Tr; the Minkowskian semiaxis < is the z line 
Imz = — tt; the point Q-^ — corresponds to z = —00; = {— (3mc)^ ~ 14.5 GeV^) corresponds to z = 0; see 
Fig. [1] If we denote a{Q'^) = F{z), RGE ([tJ can be rewritten 



^ For 13(a) = /3(a, MS) we used Pade [2/3](a) based on the known MS Cj-coefficients: C2 and C3. Using truncated (polynomial) series up 
to —j3QC-),a^ instead, changes the results almost insignificantly, by less than 1 per mil. For the quark mass values we use; rric = 1.27 
GeV and = 4.20 GeV (cf. Ref. [39]). 
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cut 



Q —plane 



z<0(0<Q<(i?) 



z— plane 



■ +1 % 



(Q=|if, ) 



(cut) 



FIG. 1: (a) Complex plane; (b) complex z plane where z = ln{Q^//i?j); the physical stripe is — tt < Imz < +7r. 



dF{z) 
dz 



(15) 



in the semiopen stripe — tt < Imz < +7r. The analyticity requirement for a(Q^) now means analyticity of F{z) 
(=> dF/dz = 0) in the open stripe — tt < Im(z) < +tt, and we expect (physical) singularities solely on the line 
Im(z) — — TT. Writing z = x + iy and F = u + iv, and assuming analyticity {dF/dz — 0), we can rewrite RGE ( 15 1 as 
a coupled system of partial differential equations for u{x, y) and v{x, y) 



du{x, y) 

dx 
du{x,y) 

dy 



Ke/3{u + iv) 



dv{x,y) 
dx 



Im/3(u + iv) 



-Im/3(u + iv) , = Re/3(M + iv) 

ay 



(16) 
(17) 



Thus, beta functions P{F) are analytic at = [ITEP-OPE condition and the expansion of I3{F) around F = 
[cf. Eq. Q] must reproduce the two universal parameters /3o and Ci — Pi/Po ("pQCD condition," where /3o = 9/4 
and ci = 16/9 for ny = 3), and solution F{z) ~ u{x,y) + iv{x,y) of RGEs (16)-(17l satisfies the initial condition 
F(0) = Cin where a — a^n is determined by Eq. (14). 



We implement high precision numerical integration of RGEs (16l-(17) with MATHEMATICA pT], for various 



Ansdtze of (3{F{z)) satisfying the aforementioned ITEP-OPE and pQCD conditions. Numerical analyses indicate 
that it is in general very difficult to obtain analyticity of F{z) in the entire open stripe — tt < Im(z) < +7r, equivalent 
to the analyticity of a(Q^) for all complex except S (— oo,0]. On the other hand, if we, in addition, require 
also analyticity of a{Q^) at = (<^ z = —oo), certain classes of /3(a) functions do give us F{z) with the correct 
analytic behavior. This = analyticity condition in general implies 



2\2i 



a(Q') = ao + ai(Q7A') + 0[(Q7A') 

where < aq = a(Q^ = 0) = F{z — — oo) < oo and ai ^ 0. Application of d/d\nQ^ 
that in the Taylor expansion of P{F) around F ~ ag the first coefficient is unity 



(18) 



P{F) = lx{F-a„)+0[{F-ao)^] , 



or equivalently^ 



P'{F)\p=a„^+l ■ 

We write our f5(F) Ansdtze in the form 

/3{F) = -/3,F\l-Y)f{Y)\Y^p^ao , 
with function f{Y) fulfilling the three aforementioned conditions 

f{Y) analytic at F (ITEP - OPE) , 
f{Y) = l + {l + c,ao)Y + OiY^) (pQCD) , 
ao;9o/(l) = 1 (Q^ = analyticity) . 



d/dz to Eq. (18) then implies 

(19) 

(20) 

(21) 



(22) 
(23) 
(24) 



If we assumed analyticity of a(Q^) in a special way, with ai = 
n > 2 and /3'{F)\p^ao = ri. This would imply ao/3o/(l) = n {> 2). 
then that in such a case the RGE solution F{z) has poles at Imz = 



in Eq. ([T8|, then we would have a(Q^) = ao + C KQ^/ A^)"] with 
Prom considerations in Appendix |a| [cf. Eqs. ( |A8| - |aTi] |] it follows 
±-K/n, i.e., Landau poles. 
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We always consider cq [= a((5^ = 0)] to be positive [note: a = {gs/2/TT)'^ > 0]. 

We will argue in more detail why and how this additional constraint [analyticity of a{Q^) at — 0] improves the 
analytic behavior of a{Q'^) = F{z) in the entire plane [z stripe), in the sense of avoiding Landau singularities. For 
this, it is helpful to consider some simple classes of beta functions which, on the one hand, allow for an implicit analytic 
solution z — G{F) of RGB (15) and, on the other hand, are representative because larger classes of beta functions 
can be successively approximated by them. Specifically, we consider f{Y) in Eq. (21) to be either a polynomial or a 
rational function'* 



fiY) 



l + J^rkY'^PiR/OhiY) , 

fe=i 

M AT 

(1 + ^ mfcy^-)/(i + J2 n^y') = m/N]f{Y) . 



(25) 
(26) 



k=l 



Here, the degrees (i?; M, N ) ar e in principle arbitrary, and the coefficients {rk] nik, ng) as well. Such Ansdtze apparently 
can fulfill all constraints (22|-(24). It is also intuitively clear that they can approximate large classes of other f3 
functions that fulfill the same constraints. 



Now we undertake the following procedure. Formal integration of RGE ( 15 ) leads to the solution 



z = G{F) , G{F{z)) 



(27) 



where Oin is the aforementioned initial value ai„ = a{Q^ — — F(Q) . Equation ( 27 1 represent s an implicit (inverted) 
equation for F ~ F{z) = G~^{z). In both cases, Eqs. (25) and (26), the integration in Eq. (27) can be performed 
explicitly. This is performed in Appendix [A] 

Here we quote, for orientation, the results for two simple examples of f{Y), a quadratic^ polynomial P[2/0]/ and 
a rational function P[l/1]/. 

In the case of quadratic polynomial we have 



fiY) = l + riY + r2Y^ , 



where ri = (1 + ciOq) due to the pQCD condition (23). The (positive) quantity = a{Q^ 
a function of the only free parameter r2 by the = analyticity condition ( 24 ) 



ao{r2) 



1 
2^ 



-(2 + ra) + v/(2 + r2)2+4ci//3o 



For the integration (27), we need to rewrite the polynomial (28) in a factorized form 



f{Y = l/t) = ^{t-h){t-t2) 



t2(r-2) 



-ri± 



4:ro 



(ri = 1 + Ciao(r2)) 



(28) 

0) is then obtained as 
(29) 

(30) 
(31) 



Integration (27) then gives the following implicit equation for F{z) = a{Q ): 



where 



(-1) 
/3o 



F{z) 



In 



ao/F{z) - 1 
ao/flin - 1 



(<i-l)(ti-t2) ' 



B2 = 



^2 



ao/F{z) 



ao/fli, 



{t2-l){t2-h) ■ 



(32) 



(33) 



* In the following we characterize such functions by the corresponding Pade-notations. 

A linear polynomial has at first only one free parameter ri = (1 + ciao) by the condition j23| ; however, this ag gets fixed by the = 
analyticity condition (|24|: ag ^ 0.1904. 
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In this solution we took into account that the coefBcient Bo/{PQao)) = 1/((1 — — i2)(/3oao)) in front of the first 
logarithm in Eq. (32| is simply unity by the = analyticity condition (24). The poles Zp, at which F{zp) = oo, 



are obtained from Eq. (32 1 by simply replacing 1/F{z) by zero 



In 



(-1) 



ao/oin - 1 



1 



2 

V B, In 



(34) 



It turns out that ap > (typically, oq ~ 0.1-0.2 and ain < 0.1). If, in addition, < r2 < rf/A, then Eqs. (31 1 imply 
ti,t2 < 0. Therefore, when < r2 < rf/A, all the arguments in logarithms in Eq. (34 1 are positive, except in the first 



logarithm where ln(— 1) = ±i7r and thus the only poles of F{z) in the physical stripe (— vr < Imz < tt) have 



(35) 



This implies that for < r2 < rf/A the considered singularity must lie on the timelike axis {Q^ < 0) and hence does 



not represent a Landau pole. We stress that for such a conclusion, the — analyticity condition (24) is of central 
importance, since it fixes the coefficient in front of ln(— 1) in Eq. (34) to be unity.^ We can derive from Eq. (34) the 
location of the pole in the plane at 



Ql. 



^i^expl^zpj 



-H^^ cxp 



-/^i„exp(Rezp 



1 



n 



flo/ lin — tj 



(36) 



On the other hand, if the aforementioned conditions are not fulfilled, we obtain — tt < Imzp < tt, representing a pole 
inside the physical z stripe and thus a Landau singularity. Specifically, when r2 < 0, we have ii > and ^2 < by 
Eqs. (31 1; numerically, we can check that in this case always oo/ain — ti > and, consequently the j — \ logarithm 



in Eq. (34) becomes nonreal and — tt < Imzp < tt, i.e.. Landau pole. 



To obse rve in more detail the occurrence and the shape of these singularities, we pursued the numerical solution of 
RGE (15), i.e., RGEs ( 16 )-(|17|), accounting for the initial condition at = (3mc)^ in the aforementioned way. In 
order to see the appearance of singularities of F{z) = F{x + iy) in the physical z stripe, it is convenient to inspect 
the behavior of |/3(F(z))| which should show similar singularities. The numerical results for |/3(F(z))|, in the case of 
r2 = and r2 = —2 are given in Figs.[2][a), (b), respectively. In these figures, we see clearly that the singularities are 



P[2/0];r, 



IACf=i;=:.y))| 




P[2/0];r, 



2.0 



l^t.F^=.y;)|. 




G 
FIG. 2: \f}{F{z))\ as a function oi z = x + iy for the beta-function \21\ with f{Y) having the form | |28| l with (a) r2 = 0; (b) r2 = —2. 

on the timelike edge Imz = ±7r in the case of r2 = where we have Oq = 1.901, ti « —1.338 [t2 is not present as f{Y) 
is a linear polynomial]. The pole moves inside the z stripe (Fe., become Landau singularities) in the case of r2 = —2, 
where we have gq — 0.5, ti w 0.756 and ^2 ~ —2.645. In Fig. pi a) we present the numerical results for the discontinuity 



® This also explains why it is nearly impossible to obtain an analytic a(Q^) if wc abandon the = analyticity condition i24\ 
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ln(cr/^2 



function pi(cr) = Ima{Q = —a — it) = lniF{z = x — ii:) = v{x,y = — tt) as a function of x = Re(z) — ii^yv^ / f^jj^^, 
for the case r2 = 0. In Fig. [sjb) the analogous curve for Rea{Q^ = ~a — it) — ReF(z — x ~ in) = u{x,y = — tt) is 
presented, for the same r2 — case. In Figs. [4] (a), (b), the corresponding curves for the r2 — —2 case are depicted. 



P[2/0]:/-2 = 

y = Pi(o-) = v(x, -n) 
0.7 r 



P[2/0]:r2 = 

T — u(x, -n) 





(b) 



X = log 



3.81' 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



X = log 



13.81- 



FIG. 3: (a) The discontinuity function pi(o-) = Ima(Q^ = —a — ie) = lmF(z = x — iir) = v(x,y = — tt) as a function of a; = Re(2) : 
ln((T//x?^), for the case when f{Y) has the form ( |28| with r2 = 0, i.e., hncar polynomial; (b) same as in (a), but for Rea(Q^ = —a — ie) ■ 
KeF{z = X — in) = u{x, y = —it). 



P[2/0]:/-2 = -2 

y = Pi(o-) = v(x, -71) 
0.7 r 



(a) 



0.6 

0.5 

0.4 

0.3 

0.2 [ 

0.1 



P[2/0]:/-2 = -2 

T — u(x, —tt) 
0.7 r 



- — (T - t€ 

fGeV-l 



= log 



(b) 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



2 4 b.81-; -8 -6 -4 

FIG. 4: Same as in Figs.js] but this time r2 = 



-0.1 



x^ log 

4 13.81- 



-2. 



We can try many other f{Y) functions, for example, the following set of functions involving (rescaled and translated) 
functions (e^ - 1)/Y and Y/{e^ ~ 1): 



EE 



(exp[-fci(r-ri)] -1) 



(exphA:2(r-r2)]-i; 



X /C(fci,Yi,fc2,1^2) 



(37) 



where the constant /C ensures the required normalization ]{Y = 0) = 1. In this "EE" case we have, at first, five 
real parameters: Oq = a(Q^ = 0) and four parameters for translation and rescaling (Yi, fci, 1^2, and k-2)- Two 
of the parameters, e.g., Yi and oq, are eliminated by conditions (23 1 and (24). We need < fci < ^2 to get 



physically acceptable behavior and fulfill the aforementioned two conditions. It turns out that, in general, increasing 
the value of Y\ tends to create Landau poles. We consider two typical cases: (1) y\ — 0.1; fci = 10; ^2 — 11; 
(2) yi = 1.1; fci = 6; fc2 = H- The numerical results for fi{F{z)) for two cases are presented in Figs. [Sja), (b), 
respectively. We see that the first case shows no sign of Landau poles, while the second case strongly indicates 
Landau poles. In Figs. [6] and [t] we present the behavior of the imaginary (v) and real (u) parts of the coupling 
F(z = X — in) — a{Q^ = —a — it) along the timelike axis of the plane for the aforementioned two EE cases. 
There is one interesting feature which can be seen most clearly in Figs. |3][a) and[6]ja): the discontinuity function 



pi(cr) = Ima((5 = —(J — it) is zero at negative Q -values above a "threshold" value: 



=) -CTthr < < 0. For 



the two cases cited there ("P[l/0]" which is "P[2/0]" with = 0, and EE with Yi = 0.1), we obtain xthr = -5.948 and 
—5.403, respectively, leading to the threshold masses Afthr = 195 MeV and 256 MeV, respectively. These threshold 
masses are nonzero and comparable to the low QCD scale Aqcd or pion mass, a behavior that appears physically 
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EE; r, = 0.1 fcj = 10. ita = 11. 



EE: y, = 1,1 Jt| = 6. = 1 1. 



|.P(F(x,y)| 4 





FIG. 5; (a) \(3{F(z))\ as a function of z = x + iy, where /3 has the form | |21| l with f{Y) having the EE form \37\ with the values of free 
parameters y\,k\,andk2 as indicated; (b) same as in (a), but with different values of parameters y\ and fci. 



EE: Yi = 0.1 ki = 10. = 11. 



y = 


Pi(o-) 


= v{x, -n) 






2.0 






(a) 


1.5 










[GeV^] 






1.0 








0.5 




8 -6 


-4 -2 




2 4 



EE: Fi = 0.1 ki = 10. &2 = H- 

y = u{x, -n) 

2.0 r 



X = log 



(b) 



1.5 



1.0 



0.5 



fGeV^l 



^ x = log 

4 b.81^ 



3.81- ; -8 -6 -4 -2 

FIG. 6: (a) The discontinuity function pi(a-) = Ima(Q^ = — cr — ie) = ImF{z = x — itt) = v{x,y = —n) as a function of x = Re(z) = 
ln((T//i? ), for the case when f{Y) is the exponential-related EE function \37l with yi = 0.1; ki = 10; k2 = 11; (b) same as in (a), but for 



Rea(Q 



■ ie) = ReF{z = x — i-w) = u{x, y = — tt). 



reasonable.^ This nonzero threshold behavior (see also Fig. [T]) for the discontinuity function pi(cr) appears because 
of the = analyticity requirement for a(Q^), Eq. (24). On the other hand, earlier, we saw that the condition 
Eq. (24) is practically a necessary condition to avoid the appearance of Landau poles of a{Q^). 



EE: Yi = 1.1 ki = 6 k2 = 11. EE: Yi = 1.1 ki = 6 kz = 11. 

y = Pi(cr) = v(x, -7r) y = u(x, -n) 

0.20 r 2.0 r 




FIG. 7: Same as in Figs.[6] but this time j/i = 1.1 and fci = 6. 



^ Furthermore, analytic couplings with nonzero Aftjn- have the mathematical property of being Stieltjes functions, and therefore their 
(para)diagonal Fade approximants are guaranteed, by convergence theorems, to converge to them as the Fade index increases I42| . 
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TABLE I: The relative deviation R[Q^] — (rhs/lhs— 1) for the Ihs and the rhs of dispersion relation ( 38 1 as obtained numerically, 
for various low positive (Q^ = 0., OT, 1.0 GeV^), for the aforementioned cases of the beta function. 







parameters 




R[Q^ = 0.0] 


R[Q^ = 0.1] 


RIQ"" = 1.0] 


P[2/0] 




ra = 0.0 




3.3 X 10"^ 


4.6 X lO""* 


7.0 X 10"^ 


P[2/0] 




r-2 = -2.0 




-0.62 


-0.38 


-0.09 


EE 




= 0.1, fci = 10.0,^2 


= 11.0 


4.7 X 10"^ 


4.8 X 10"^ 


6.5 X 10"^ 


EE 




= 1.1, fci = 6.0, fc2 


= 11.0 


-0.82 


-0.68 


-0.19 



While Figs. [2] and [5 provide only a visual indication of whether the coupling a{Q^) is analytic, there is a more 
quantitative, numerical test for the analyticity. Namely, application of the Cauchy theorem_ implies for an analytic 
a{Q^), with cut along the negative axis 
starts effectively at cr = ath: 



ML 



< —AlLi the well-known dispersion relation (|6|) where the integration 



da 



Art 



(38) 



where pi(cr) = lma{Q^ = —a — ie). The high precision numerical solution of RGE (15) gives us a{Q^) — F{z) in the 
entire complex plane, incl uding the negative semiaxis. This allows us to compare numerical values of the Ihs and 
rhs of dispersion relation (38), for various values of . 



It turns out that, for low positive < 1 GeV^, the numerical uncertainties of the obtained results for the rhs 
of Eq. ([38| are of the order of per cent (using 64-bit MATHEMATICA [41] for Linux), and they slowly increase 
with increasing . If the deviation of the rhs from the Ihs is more than a few percent, then this represents a strong 
indication that the resulting a{Q'^) is not analytic. In Tablejljwe present the relative deviations for the aforementioned 
two P[2/0] and the two EE cases. Inspecting these deviations, we can clearly see that a{Q^) in the P[2/0] case with 
r2 = —2 and the EE case with Yi = 1.1 is nonanalytic; in the other two cases, the table gives strong indication that 
a{Q^) is analytic. 



III. EVALUATION OF LOW-ENERGY OBSERVABLES 



The semihadronic r decay ratio is the most precisely measured low-energy QCD quantity to date. The measured 
value of the "QCD-canonical" part r-r — a + ©(a^), with the strangeness and quark mass effects subtracted, 
„(™p) 



0.203 ±0.004 (cf. Appendix |B[). Experimental values of other low-energy observables, such as (spacelike) sum 
rules, among them the Bjorken polarized sum rule (BjPSR) dBj(Q^), are known with far less precision. The minimal 
analytic (MA) model [IHHl [20], with the value of A such that high-energy QCD observables are reproduced, turns out 
to give for this quantity too low values ~ 0.14 [51 [55] unless the (current) masses of the light quarks are taken to 
be unrealistically large {iriq « 0.25-0.45 GeV) or strong threshold effects are introduced [33]. Further, MA does not 
fulfiU the ITEP-OPE condition (g since ja^'^^HQ^) _ apt{Q^)\ - (A^/Q^). 

The approach described in the previous Sec. [llj automatically fulfills the ITEP-OPE condition ([5]); however, the 
analyticity of a(Q^), i-e., the absence of Landau poles, is achieved only for limited regions of the otherwise free 
parameters of the /3 function. For general anQCD models, the evaluation of massless spacelike observables 'D{Q^) 
such as BjPSR and Adler function, and for the timelike observable r,-, is presented in the sequence of Appendixes [C| 
|D][E) particularly Eqs. (|E9|)-( |eT2| ) for spacelike and ( [E22| )-( |E25| ) for r^. In the cases considered in this work, the beta 
function /3(a) is analytic at a = (due to the ITEP -OPE condition), and therefore the higher order analogs An+i in 
those Appendixes are simply An+i = a"-'^^, cf. Eq. (C29). Furthermore, here we use all t he t ime the notation Ai = a 
for the analytic coupling, and An+i = On+i ioi the logarithmic derivatives of a [cf. Eq. (C5|]. 

In Table [TTj we present the resulting values of RSch parameters C2, C3 and C4 [cf. Eq. ([4|j, for some typical choices of 
input parameters in four forms of f{Y): P[l/0], P[3/0], P[l/1], and EE. Here, P[M /N] is the general notation for Pade 
form Eq. (Al) in Appendix [A{ P[Af/0] is thus a polynomial of degree M; EE is the Ansatz (37l involving exponential 
functions. The otherwise free parameters ("input") of the models are chosen such that the analyticity is maintained, 
i.e., no Landau poles. The case P[l/0] is in fact the aforementioned case of P[2/0] with r2 = 0, cf. Eq. (28), and it 
has no free parameters. The cases P[3/0] and P[l/1] have each one free input parameter; for P[3/0] the first root ti 
is the specified input, and for P[l/1] the first pole wi, where the notation (Al) of Appendix |A] is used. The case EE 
is given in Eq. (37), and has three free par ame ters. We recall that an apparently additional parameter in the Ansatze 
for f{Y) is fixed by the pQCD condition (23). In addition, we present the values of a{Q^) at the initial condition 
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TABLE II: Four cases of /3 function {f{Y)), with chosen input parameters. Given are the resulting RSch parameters c„ 
(n = 2,3, and 4), and the values of a{Q^) at — (3mc)^ and = 0. Further, the resulting threshold parameter Xthi and 
the threshold mass Afthr (in GeV) are given. Recall that a((3mc)^MS) = 0.07245. 



/ 


Input 


C2 


C3 


C4 


a((3mc)^) 


ao = a(0) 


Xthv 


Mth. [GcV] 


P[l/0] 




-37.02 








0.06047 


0.1901 


-5.948 


0.195 


P[3/0] 


ti = 1 + iO.45 


-39.55 


115.88 


-105.80 


0.06066 


0.4562 


-11.092 


0.015 


P[l/1] 


Ml = "0.1 


-37.54 


18.84 


-9.46 


0.06048 


0.1992 


-6.060 


0.184 


EE 


Yi = 0.1, fci = 10.0, fcz = 11.0 


-10.80 


-157.62 


-644.32 


0.06544 


0.2360 


-5.403 


0.256 



scale fif^ — (3mc)^ {rric = 1.27 GeV) and at = 0; and the threshold value a;thr of the discontinuity function 
Pi{a) — lma{—a — ie), where: zthr — Xthi — i^^, fthr — (3mc)^ exp(a:thr)- Further, the corresponding threshold mass 
Mthr is given [Mthr = Sm^ exp(a;thr/2)]. 

For two of these models (P[l/0], and EE), we depict in Figs. |8] and [o] the form of f{Y) and j3{x) functions for real 
values oi Y = a/aQ and positive values of x = a > 0, respectively. In Figs. lOpT we present the running coupling 
a{Q^) as a function of for positive in the two models; there we include, in addition, the higher order analytic 
couplings a-n+i (n = 1, 2). 




FIG. 8: (a) f{Y) function as defined by Eq. | [2l) , for real values of y = a/a(0), for the case of / being P[l/0] linear function (<t4> P[2/0] 
with r2 = 0); (b) function for the same case, for positive x = a. 




The model with / = P[l/0] is, at first sight, very similar to the model of Ref. ^43j which was obtained on the basis 
of the principle of minimal (renormalization scheme) sensitivity (PMS) [33] applied to the QCD part of Rf,+e~{s) 
ratio. There, the beta function is also a polynomial of the fourth degree, i.e., f{Y) is linear, and it ha s a fi nite positive 
value of a{Q^ = 0) = Cq. It turns out that for the beta function of Ref. |43j the conditions (22) and (23) are fulfilled, 
but not the condition oi Q-^ — analyticity Eq. (24). As argued in the present paper, such beta function will give 
unphysical (Landau) poles, although in this case not on the positive axis. Specifically, for n/ = 2 and n/ = 3 the 
= analyticity condition (24) yields in the P[l/0] case the values ao — 0.1761 and 0.1901, respectively, while the 
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-0.05 
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AQ 
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14 



0.20 



0.15 



0.10 



0.05 



0.00 



-0.05 ^ 
0.0 



P[l/0] 

m 



(b) 



5aM\ . 



25ji3(e:). 



0.2 



0.4 



0.6 
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1.0 



Q- [GeV^] Q- [GeV-] 

FIG. 10: (a) Analytic coupling a(Q^) and its higher order analogs a,i_|_i (n = 1,2) as defined in Eq. jC5| , for positive Q^, for the model 
P[l/0]. For better visibility, the higher order analogs are scaled by factors of 5 and 5^, respectively, (b) Same as in (a), but at lower Q^. 
We recall that, formally: a„+i = a"+i + 0{a"+^). 
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EE: Yi = 0.1 ki = 10. kj = 11. 
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FIG. 11: Same as in Fig.lTol but for the model EE, Eq. (IstI 



0.4 0.6 
[GeV^] 



0.8 
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values of ao in Ref. [331 are ap = 0.263 and 0.244, respectively. We checked numerically that this PMS solution leads 
to (Landau) poles of a{Q^) at w (-0.027 ± i0.065) GeV^ for nj = 2, and at w (-0.031 ± i0.032) GeV^ for 
n/ = 3 (massless quarks). 

Let us now apply these results to calculating low-energy QCD observables. 

We start with r-r- 



In Table III we present the predicted values of rr for the choices of /? functions and input parameters given in Table 
[ll| Therein we separately give (in each line) the four terms of the truncated analytic series for and then quote their 
sum. Furthermore, for each model of f{Y) we present the results for basically two different ways of treating the higher 
orders. In the first row of each model, the results of the series (E22| are presented, which performs leading- /3o (LB) 



resummation and adds the (three) beyond-the-leading-/3o (bLB) terms organized in contour integrals of logarithmic 
derivatives a„+i {n = 1,2,3). In the second line, the an alogo us results are presented, where now the (three) bLB 
terms are contour integrals of powers An+i = a"+^, Eq. (E24). At each of the entries, the corresponding terms are 
given when no LB resummation is performed, cf. Eqs. (E23l, (E25). The RScl parameter used is C = 0, i.e., the 
radius of the contour in the plane is m^. In the last column, the relative variation of the sum is given when the 
RScl parameter is increased from C = to In 2, i.e., the radius of the contour integration is increased to 2m^. The 
results using the powers a""*"^ for the bLB (or: higher order) contributions show significantly less stability under the 
RScl variation; the reason for this lies in two numerical facts: 

• The expansion coefficient (iAdOs of the latter series is usually larger than the corresponding coefficient (TAdOs of 
the series containing a„+i: KtAdOal > KTAdOal; this seems to be true in all the RSch's dictated by the presented 
/3 functions. 



Apparently in all cases we have |a„+i| < 



I, although formally a„+i = a"+i -I- C'(a"+2). 
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TABLE III: The four terms in truncated analytic expansions (E22l and (E24l for rw, i.e 
the three bLB terms organized in contour integrals of An+i = 5; 

In pa renth eses are th e corresponding results when no LB resummation is performed, i.e 
Eqs. (E23l and (E25l, respectively. The RScl parameter is C = 0. The last column contains variations of these truncated sums 



when the 



with LB contributions resummed and 
1 (first line) and of ^„+i — a"^^ (second line of each model). 

the truncated analytic expansions 



RScl parameter C increases from to In 2. 



/ 


: LB (LO) 


NLB (NLO) 


N^LB (N^LO) 


N^'LB (N=^LO) 


Sum (sum) 


S (C dependence) 


P[l/0] 


0.1135(0.0940) 
0.1135(0.0940) 


0.0006(0.0123) 
0.0007(0.0137) 


0.0139(0.0214) 
0.0209(0.0340) 


0.0007(0.0012) 
0.0091(0.0113) 


0.1287(0.1289) 
0.1442(0.1529) 


-0.2%( 
-2.8%( 


-0.4%) 
-2.7%) 


P[3/0] 


0.1200(0.0954) 
0.1200(0.0954) 


0.0007(0.0131) 
0.0007(0.0141) 


0.0184(0.0275) 
0.0233(0.0369) 


-0.0009(0.0000) 
0.0067(0.0087) 


0.1381(0.1360) 
0.1507(0.1550) 


-0.3%( 
-2.4%( 


-0.8%) 
-2.9%) 


P[l/1] 


0.1142(0.0941) 
0.1142(0.0941) 


0.0006(0.0124) 
0.0007(0.0138) 


0.0146(0.0224) 
0.0213(0.0344) 


0.0005(0.0011) 
0.0088(0.0109) 


0.1300(0.1300) 
0.1450(0.1532) 


-0.2%( 
-2.8%( 


-0.5%) 
-2.7%) 


EE 


0.1348(0.1088) 
0.1348(0.1088) 


0.0009(0.0173) 
0.0009(0.0180) 


0.0025(0.0156) 
0.0033(0.0224) 


0.0048(0.0061) 
0.0102(0.0173) 


0.1466(0.1478) 
0.1528(0.1666) 


-0.8%( 
-2.8%( 


-1.2%) 
-3.7%) 



jola rized sum rule (BjPSR) results dBj(Q^) for the four considered /3 Ansdtze, evaluated with the truncated 
with LB contributions resummed and the three bLB terms oc An+i = 5„+i (first line) 



TABLE IV: Bjorken 
analytic expansions (E9l and (Elll, i.e 

and oc An+i = a""*"^ (second line) . In parentheses are the corresponding results when no LB resummation is performed, i.e 



truncated analytic expansions Eqs. (ElOl and (E12l, respectively. The RScl parameter is C = 0. In brackets, the corresponding 
variations of the results under the RScl variation are given (see the text for details). For explanation of the experimental values 
in the last (four) lines, see the text for details. 



/ 


dBj(Q') : 


= 1.01 GeV^ 


= 2.05 GeV' 




0^ = 2.92 GeV^ 




P[l/0] 


0.1343[+0.3%] 
0.1535[-4.1%] 


(0.1420[ 
(0.1974[ 


-1.9%]) 
-5.1%]) 


0.1208[-0.1%] (0.1255[- 
0.1313[±2.8%] (0.1552[- 


-0.5%]) 
-4.1%]) 


0.1140[-0.2%] (0.1173[- 
0.1218[-2.4%] (0.1393[- 


-0.7%]) 
-3.6%]) 


P[3/0] 


0.1609[-0.4%] 
0.1773[-3.7%] 


(0.1630[ 
(0.2053[ 


-1.9%]) 
-6.3%]) 


0.1366[-0.4%] (0.1361[- 
0.1456[-2.6%] (0.1587[- 


-2.0%]) 
-4.6%]) 


0.1261[-0.4%] (0.1249[- 
0.1329[-2.2%] (0.1417[- 


-1.9%]) 
-4.0%]) 


Pll 


0.1373[+0.2%] 
0.1561[-4.0%] 


(0.1450[ 
(0.1985[ 


-1.5%]) 
-5.3%]) 


0.1226[-0.1%] (0.1270[- 
0.1329[-2.8%] (0.1557[- 


-0.6%]) 
-4.2%]) 


0.1154[-0.2%] (0.1184[- 
0.1231[-2.4%] (0.1396[- 


-0.9%]) 
-3.7%]) 


EE 


0.1507[+0.3%] 
0.1436[+0.7%] 


(0.1659[ 
(0.2300[ 


-3.7%]) 
-6.8%]) 


0.1338[±0.1%] (0.1434[- 
0.1304[±0.5%] (0.1725[- 


-1.0%]) 
-5.1%]) 


0.1256[+0.1%] (0.1324[- 
0.1232[+0.4%] (0.1521[- 


-1.0%]) 
-4.4%]) 


Exp. (a): 


0.23 ±0.18 






0.11 ±0.11 




0.09 ±0.07 




= -0.040 ± 0.028 


0.23 ±0.12 ±0.13 




0.11 ±0.09 ±0.06 




0.09 ±0.05 ±0.05 




Exp. (b): 


0.30 ±0.18 






0.15 ±0.11 




0.11 ±0.07 




= -0.024 ± 0.028 


0.30 ±0.12 ±0.13 




0.15 ±0.09 ±0.06 




0.11 ±0.05 ±0.05 





Furthermore, the variations of the result under variations of RScl are generally smaller when LB resummation is 
performed. Therefore, we will consider as our preferred choice the evaluated values of the first lines (not in parentheses) 

the evaluations using a„+i for the higher order contributions, i.e., Eq. (E22). 
"sum" 



of each model in Table III 



We note that the obtained values of (see the 



experimental value r'f^^'^ 



in Table III I are all much too low when compared with the 



0.203± 0.004 (cf. Appendix |B[). In fact, the free parameters in the Ansdtze for f{Y) of the 



beta function were chosen in Tables IIpII in such a way as to (approximately) maximize the result for Tt while still 
maintaining analyticity of a{Q^) (i.e., no Landau singularities).^ We can see that the preferred evaluation method, 
i.e., the first line of each case, gives us always a value < 0.15. We tried many choices for the function f{Y) of 
Eq. (21), fulfilling all conditions (22)-(24|, and scanning over the remaining free parameters in fiY). 
that rr < 0.16 always as long as Landau poles were absent.^ 
poles appeared, was it possible to increase r^. beyond 0.16. 



It turned out 

Only when free parameters were chosen such that Landau 



As the s econ d example we consider the Bjorken polarized sum rule (BjPSR) dBjiQ )■ 



In Table 



IV we present results for c?Bj(Q^) in the aforementioned cases, at three of those low values of where 



When f{Y) is P[2/0], it turns out that the largest evaluated value of r-r is obtained when r2 = in Eq. ( |28| , i.e., when f{Y) reduces 
to a linear function P[l/0]. 

In some cases, e.g., when increasing the value of Yi in the case EE, the preferred evaluation method, Eq. ( |E22[ |, gives us values of 
between 0.15 and 0.16. However, in such cases, it is not any more clear that the analyticity is maintained; increasing Yi even further 
leads to clear appearance of Landau poles. 
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experimental results are available: = 1.01, 2.05, and 2.92 GeV^. As in the previous Table III the first line of 



each model contai ns th e results with our preferred method, i.e., LB resummation and usage of an+i for the bLB 
contributions, Eq. (E9); the second line represents the results of LB resummation and the usage of a""*"^ powers for 



the bLB contributions, Eq. (Ell). In the parentheses, the corresponding results are given when no LB resummation 
is performed, Eqs. (ElO) and (E12), respectively. In the corresponding brackets, the variations of the results are given 
when the RScl parameter varies either from C = (/i^ = Q^) to C = ln2 {p^ = 2(5^), or from C = to C = ln(l/2) 
(/i^ = Q^/2) - the larger of the variations is given. As in the case of r^., we see that the most stable evaluation under 



variations of RScl is the LB resummation and the usage of a„+i for the bLB contributions, Eq. (E9) 



For comparison, we include in Table pV| (last lines) three sets of experimental data based on the JLab CLAS EGlb 
(2006) measurements [45| of the T\~'^(Q'^) sum rule for spin-dependent proton and neutron structure functions g^''" 
[46] . r^~" is connected to dej in the following way: 



1 



(Q' 



9A_ 

6 



dxBj (.gi(a;Bj,Q^) -5?(a;Bj,Q^)) 

(i-rfB.(Q^))+E%y&^^ 



(39) 
(40) 



J=2 



where qa = 1.267 ± 0.004 [39j is the triplet axial charge, 1 — c?Bj(Q'^) ~ ^ — a(Q^) + C'(a^) is the nonsinglet leading- 
twist Wilson coefficient, and IJ^^^ /Q'^-'~^ {j > 2) are the higher-twist contributions. If we take into account the 
data with the elastic contribution excluded, we can restrict ourselves to the first higher-twist term /i4~"/(5^. The 
elastic contribution affects largely only the other higher-twist terms ~ 1/(Q^)''~^ with j > 3, as has been noted in 
Refs. [47l |48]. Moreover, the exclusion of the elastic contribution leads to strongly suppressed higher-twist terms 
~ 1/(Q2)^"^ with j > 3 [17] in pQCD and MA (APT) approaches. The first experimental set (a) for dBjiQ'^) in Table 



IV is obtained from the measured values of r^~"(Q^) (with the elastic part excluded) by subtracting the "/Q^ 
contribution as obtained by a 3-parameter pQCD fit [45]: /^4~" ~ M4~'^('9 = IGeV) = — 0.040±0.028;^° the second set 
(b) is obtained in the same way, but now by subtracting the /i4~"/Q^ contribution obtained by a 4-parameter pQCD 
fit @S]: /X4~'^ « /^4~"(Q = IGeV) = — 0.024 ± 0.028. In the second fine of each experimental set, the uncertainties 
were split into the contribution coming from the uncertainty of the measured value of T^~'^{Q^) and the one from the 
uncertainty of the fitted value /i4~" [45 . 

We see from Table [TV] that the evaluated values for BjPSR lie in general relatively close to the central experimental 



values (iBj(Q )exp.: dsjiQ )c 



0.23 (or 0.30) for = 1.01 GeV^ 0.11 (or 0.15) for = 2.05 GeV^ 0.09 (or 



0.11) for = 2.92 GeV . However, in contrast to r,-, the experimental uncertainties are now much larger and the 
theoretical predictions lie well within the large intervals of experimental uncertainties. 



IV. TACKLING THE PROBLEM OF TOO LOW 



The problem of too low r^., encountered in the previous Section, appears to be common to all or most of the anQCD 
models. For example, in the MA of Shirkov, Solovtsov and Milton [3H3 [23111]) when adjusting A to such a value as to 
reproduce higher energy QCD observables (Q^ > 10^ GeV^), i.e., A w 0.4 GeV, the resulting^^ value of (massless and 
strangeless) is about 0.140-0.141 [SKTHj^S], much too low. The results of the previous section indicate that this 
problem persists even in anQCD models which, unlike MA, fulfill the ITEP-OPE condition The aspect of anQCD 
models which appears to cause the tendency toward too low values of Vr is the absence of (unphysical) Landau cut 
along the positive axis (0 < < Qlp)-^^ Therefore, we are apparently facing a strange situation: 



" Almost the same value was obtained by the authors of Refs. |43|48]: /il "/M^ ^ -0.048 corresponding to fi^ " -0.042 (Ref. 1470 . 
and /^J^~"/A/p —0.042 corresponding to /^4~" ~ —0.037 (Ref. 1481 , accounting for the Q^-dependence of /^4~° due to RG evolution.). 
The interesting aspect is that they applied MA (i.e., APT) model of Refs. [1] [5] in the fit of the aforementioned JLab data, then 
obtaining the 1/Q^-term as the sum of the contribution from the MA (APT) series and the contribution of the explicit l/Q'^-term 
(obtained through fit). Such a sum of 1/Q^-terms, in their model, is not interpreted as originating entirely from the IR regime since 
MA does not satisfy the conditions of Eq. 

The value A = 0.4 GeV corresponds to the A value in the Lambert function 1491 for the (MA) coupling ^i(Q^) in the 't Hooft RSch 
Axjambort = 0.551 GeV. In general, it can be checked that the following relation holds: A^ambert ~ Aexp(0.3205), and this holds 
irrespective of whether we consider pQCD or MA couplings. 

A somewhat similar reasoning can be found in Ref. [30]. 
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TABLE V: Four models of /3 function (/(V)) of the previous section, with modification Eqs. ( |41[ )-(42 1, with inputs as given 
in Table [ll] and the values of the additional input parameters K and B {1 ^ K ^ B) adjusted so that the evaluation 
. _ . 0.203. Given are the resulting RSch parameters c„ {n = 2, 3, and 4), and the values of a{Q^) at 

and the threshold mass Mthr (in GeV). 



method Eq. ( E22 1 gives r-r 

= (3mc)^ and — 0, as well as the resulting threshold parameter xth 



/old 


Input /fact 




C2 


C3 


C4 


a((3m,)2) 


ao = a(0) 


Xthr 


A/thr [GeV] 


P[l/0] 


B = 4000, K = 


6.71 


-222.06 


-329.13 


2.047 X 10^ 


0.05763 


0.1904 


-6.331 


0.161 


P[3/0] 


B = 5000, K = 


44.5 


-249.65 


-260.93 


5.036 X 10** 


0.05430 


0.4597 


-12.023 


0.009 


P[l/1] 


B = 4000, K = 


7.11 


-216.04 


-298.77 


1.799 X 10^ 


0.05761 


0.1995 


-6.448 


0.152 


EE 


B = 1000, 7^ = 


5.4 


-106.80 


-326.71 


1.721 X 10*^ 


0.06125 


0.2370 


-5.887 


0.201 



• In pQCD the Landau cut of the coupUng gives a numerically positive contribution to r,-, and pQCD is able to 
reproduce the experimental value of (cf. Refs. [301 EDHS], because of this (unphysical) feature of the theory. 

• In anQCD the physically unacceptable low-energy (Landau) singularities of the coupling are eliminated, but 
then the values of tend to decrease too much. 



Here we indicate one possible solution to this problem (cf. also our shorter version [31]). Table [Till indicates that 
the LB-resummed contribution to r,- cannot surpass the values 0.14-0.15. We performed many trials with various 
forms of f{Y) functions and were not able to obtain larger values of r'^^\ But the N^LB term, which is the only 
nonnegligible bLB term in Table III can be increased by increasing the coefficient (TAdi)2 of expansion ( |E22[ ) while 
maintaining, at least approximately, the values of a{Q^) and Un+iiQ^) for most of the complex Q^. It can be deduced 
from the presentation in Appendix [e| that the RSch dependence of coefficient (TAdi)2 is in the contribution (— C2 + C2). 
Therefore, if we multiply the f(Y) function by a factor /fact(i^), which is close to unity for most of the values of Y 
(= a/oo) but which significantly decreases the RSch parameter C2, the value of (rAdi)2 will increase while the values 
of of a[Q'^) and a„+i((5^) will not change strongly for most of the complex values. This can be achieved by the 
following replacement: 



/old(r) /„ew(r) - 

with : /fact(i^) = 



/old(F)/fact(l') : 

[l + BY^) 
{1 + {B + K)Y^ 



(1< X < B) . 



(41) 
(42) 



The function /fact(^) is really close to unity for most y's because K <^ B\ and it decreases the C2 RSch parameter 
to low negative values [cf. Eq. (|4|] because 1 <C if (c2 ~ More specifically, expansion in powers of y = a/ao 

then gives the RSch coefficients c„ with large absolute values C2 ~ —K/a'^{^ C3 w —ciK/aQ{^ C4 ~ 

BK/a^{^ BK)\ etc. This i mplies that the coefficients (TAdi)n, (^AdOn, (^AdOn and (dAdOn appearing in analytic 
expansions Eqs. (E20 )-( [E25 ) behave as ~ — C2 ^ K for n = 2; ~ ±02,-03 ~ iif for n = 3; — C4 —BK for 
n = 4; etc. Therefore, these coefficients are large for n = 2, 3, and even much larger for n > A. In fact, it turns out 
that the larger B is, the less the LB contribution r)- decreases. However, then the absolute values of coefficients 
of analytic expansions Eqs. (E20)-(E25) incr ease e xplosively for n > 4. On the other hand, when B {'3> 1) decreases, 
the aforementioned divergence of the series (E20) at n > 4 becomes less dramatic, but then decreases and it 

becomes difficult to reproduce the experimental value 0.203. We chose the values of B in each model such that, 
roughly, rj- « 0.10 or above (if possible). 

Further, it turns out that these modifications (i.e., inclusion of /tact) do not destroy the analyticity of a{Q'^). The 
(two- and three-dimensional) diagrams presented in the figures of the previous section change only little when the 



modification factor (42) is introduced in the corresponding beta functions. 

The numerical results in the models of Tables|ll III IV of the previous section, modified by replacements (41 1-(42) 
in the aforementioned way so that the preferred evaluation method Eq. (E22| gives — 0.203, are given in the 
corresponding Tables [V VI |VII 



When comparing Table VI with Table HI we see that the modification ( 41 )-( 42 ) really results in a significantly larger 



N LB contribution (and a somewhat larger N LB contribution) to r-r, reaching in this way the middle experimental 
value Tt = 0.203. The variations S under the variations of RScl are now larger in Table [Vl| than in |III| nonetheless. 



The next-to- leading- /3o (NLB) term cannot be increased in this way, because the coefficient (TAdi)i = 1/12 turns out to be RSch 
independent (and small). 
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TABLE VI: The evaluated quantity as in Table III but now with modifications Eqs. (41l-(42|, as given in Table [v} so that 
the evaluation method Eq. (E22l gives — 0.203. 



/old — ///fact 


: LB (LO) 


NLB (NLO) 


N^LB (N^LO) 


N^LB (N^LO) 


Sum (sum) 


5 (C dependence) 


P[l/0] 


0.1060(0.0880) 
0.1060(0.0880) 


0.0006(0.0110) 
0.0006(0.0121) 


0.0907(0.0974) 
0.1264(0.1373) 


0.0057(0.0063) 
0.0552(0.0438) 


0.2030(0.2026) 
0.2882(0.2812) 


-1.4%(^ 
-8.4%(- 


-1.5%) 
-10.1%) 


P[3/0] 


0.0997(0.0815) 
0.0997(0.0815) 


0.0005(0.0099) 
0.0005(0.0104) 


0.0967(0.1029) 
0.1143(0.1230) 


0.0061(0.0068) 
0.0447(0.0347) 


0.2030(0.2011) 
0.2592(0.2496) 


-2.5%(- 
-7.6%(- 


-2.7%) 
-9.6%) 


P[l/1] 


0.1064(0.0880) 
0.1064(0.0880) 


0.0006(0.0111) 
0.0006(0.0121) 


0.0902(0.0971) 
0.1229(0.1338) 


0.0058(0.0063) 
0.0532(0.0423) 


0.2030(0.2025) 
0.2832(0.2762) 


-1.6%(- 
-8.3%(- 


-1.7%) 
-10.0%) 


EE 


0.1247(0.0987) 
0.1247(0.0987) 


0.0007(0.0146) 
0.0008(0.0149) 


0.0678(0.0786) 
0.0787(0.0934) 


0.0097(0.0108) 
0.0432(0.0385) 


0.2030(0.2027) 
0.2474(0.2456) 


-2.4%(- 
-8.8%(- 


-2.8%) 
-10.3%) 



TABLE VII: The evaluated quantity BjPSR dsjiQ^) as in Table IV but now with modifications Eqs. (41h(42l, as given in 
Table |V] The experimentally measured values are given in the last four lines of Table IV (see the text there for details) . 



/old = //.ffact 


dBj(Q^) : = 1.01 GeV^ 


Q2 = 2.05 GeV^ 


Q2 = 2.92 GeV^ 


P[l/0] 


0.2138[-2.9%] (0.2199[-3.5%]) 
0.3795[+15.0%] (0.3673[+15.2%]) 


0.1895[-1.7%] (0.1927[-2.0%]) 
0.2803[+12.4%] (0.2697[+12.9%]) 


0.1761[-2.2%] (0.1782[-2.6%]) 
0.2442[+ll.l%] (0.2344[+11.8%]) 


P[3/0] 


0.2485[-4.9%] (0.2476[-5.7%]) 
0.3485[+14.8%] (0.3221[+16.0%]) 


0.2008[-4.5%] (0.1991h5.3%]) 
0.2579[+11.5%] (0.2392[+12.8%]) 


0.1813[-4.3%] (0.1795[-5.0%]) 
0.2252[+10.2%] (0.2093[+11.5%]) 


Pll 


0.2185[-2.1%] (0.2244[-2.5%]) 
0.3742[+15.2%] (0.3618[+15.4%]) 


0.1909[-2.1%] (0.1938[--2.5%]) 
0.2761[+12.3%] (0.2654[+13.0%]) 


0.1767[-2.5%] (0.1785h3.0%]) 
0.2406[+11.0%] (0.2308[+11.8%]) 


EE 


0.2166[-3.0%] (0.2281[-4.1%]) 
0.3246[+18.4%] (0.3416[+18.3%]) 


0.1879[-2.3%] (0.1938[-3.0%]) 
0.2380[+13.6%] (0.2421[+14.6%]) 


0.1728[-2.7%] (0.1765[-3.6%]) 
0.2074[+11.6%] (0.2081[+12.8%]) 



the evaluation method of Eq. (E22| is still the most stable under the RScl variations. However, now the series for 
Tr is strongly divergent when terms N^LB and higher are included, for th e rea sons m entioned earlier in this section. 
For example, the N^LB contribution to Vt, in the methods of Eqs. (E22| and (E23l which use a„+i in higher order 
contributions, is estimated to be ^ —10" ~ —1. Specifically, when the RScl parameter is C = 0, these terms are 
estimated to be -3.1 (P[l/0]); -2.0 (P[3/0]); -3.7 (P[l/1]); -1.0 {EE)}^ 

It remains unclear how to deal with such an analytic series, which has relatively reasonable convergence behavior 
in its first four contributions and behaves uncontrollably for n > 4. One might consider this behavior as an indication 
of the asymptotic series nature of the expansion ("precocious asymptoticity" ) . Certainly, this divergence problem 
appears to be the price that is paid to achieve in anQCD the correct value r,- sa 0.20 via /3 function modification 
Eqs. (41|-(42). The modified beta functions /3(a) now acquire poles and zeros on the imaginary axis close to the 
origin in the complex a plane: Opoio = iia{0) / y/ B + K, a^cm = ±m(0)/v'S. Consequently, the convergence radius 
of the perturbation expansion of /3(a) in powers of a becomes short: R = a{0)/\/B + K. Nonetheless, /3(a) remains 
an analytic function of a at a = 0, fulfilling thus the ITEP-OPE condition ([5]). We note that such a modification 
of the beta function brings us into an RSch where the absolute values of the (perturbative) RSch parameters c„ rise 
fast when n increases. There is no physical equivalence of such RSch's with the usual RSch's such as MS or 't Hooft 
RSch (where c„ = for n > 2). For example, in these two latter RSch's, the coupling a(Q^) is not even analytic. 
Physical nonequivalence can even be discerned between, on the one hand, the much "tamer" RSch's of the previous 
Section which give analytic a(Q^) (see Table [ll| and, on the other hand, the aforementioned nonanalytic RSch's MS 
or 't Hooft. 

When comparing the evaluated BjPSR values for the beta functions modified by Eqs. (41)-(42|, as presented in 



Table |VIH with those of unmodified beta functions as presented in Table |IV[ we note that the modification increases 
the values of BjPSR, generally to above the experimental middle values. Nonetheless, the results generally remain 
inside the large intervals of experimental uncertainties. The variations of the results under the variation of the RScl 



^* When using evaluation methods of Eqs. | |E24| | and | |E25| which use powers a'^+-'^ instead, these estimated terms are: -22.9 (P[l/0]); 
-3.9 (P[3/0]);-20.1 (P[l/1]); -2.9 (EE). These terms have significantly higher absolute values than those for the methods of Eqs. ( |E22| | 
and (E23l, although the estimated coefficients are the same in both cases. The reason for this difference lies in the fact that |a^(Q'^)| > 
1^5 (Q^) I for most values of (complex) Q-^. It appears to be a general numerical fact in all models presented in this work that |a"+^ (Q^)l > 



|<in+l(Q )| > 1)> although formally a^+i 



+ Cl(a"+^). 
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are now larger. 

The evaluation methods of Eqs. ( E9 1 and ( ElO I , for spacelike observables such as BjPSR, and the analogous methods 
of Eqs. (E22) and (E23) for the timelike r^., whic h use lo garithmic de rivati ves Qn +i; are significantly more stable under 
the variation of RScl than the methods of Eqs. (Ell), (E12|, (E24| and (E25I, which use powers a"+^. This can be 



seen clearly by comparing the variations (percentages) of the first and the second line of each anQCD model in Tables 
VI and VII In this sense, the method of Eqs. (E9) for spacelike, and ( E22[ ) for timelike observables, which performs 
LB resummation and uses logarithmic derivatives a„+i for the bLB contributions, remains the preferred method, as 
in the previous section. 

We wish to add a minor numerical observation. Unlike the results of the previous section where the LB resummation 
improved significantly the stability under the RScl variation, this improvement becomes less clear in the results of 
the present section, as can be seen by comparing the variations (percentages) outside the parentheses with the 
corresponding one s in s ide the parentheses. This can be understood in the following way: the modification of /3 
functions by Eqs. (41|-(42| introduced, via large values of |c„|'s, in the expansion coefficients (i„+i and of the 

(spacelike) observables (here the Adler function and BjPSR) numerically large contributions w ~Cn+i/n which are 
not a large-/3o (LB) part of these coefficients. The latter is true because the LB part of d„+i and dn+i is ^ /^o while 
Cn+i = Pn+i/Po ~ /3o ('^f- Appendixes [P] and [e|) . Therefore, the LB parts of the coefficients are now not dominant, 
and the LB resummation cannot be expected to improve significantly the RScl stability of the result. 



CONCLUSIONS 



In this work we tried to address two aspects which are not addressed by most of the analytic QCD (anQCD) models 
presented up to now in the literature: 

• Several anQCD models, in particular the most widely used anQCD model (minimal analytic: MA) of Shirkov, 
Solovtsov, and Milton [31-1^ [20]. give significantly too low values of the well-measured (QCD-canonical) semi- 
hadronic r-decay ratio once the free parameter(s) (such as A) are adjusted so that the models reproduce the 
experimental values of high-energy QCD observables {\Q'^\ > 10^ GeV^), cf. Refs. [5t 128). 

• In most of the anQCD models presented up to now, the ITEP-OPE condition ([s]) is not fulfilled. Hence 

such models give nonperturbative power contributions ~ (A /Q"^)^ of ultraviolet origin in the (leading- twist 
part of the) spacelike observables 2?(Q^), contravening the ITEP-OPE philosophy [231 Hi] which postulates that 
nonperturbative contributions have exclusively infrared origin. If the latter philosophy is not respected by a 
model, application of the OPE evaluation method in such a model becomes questionable. 

In this work, the second aspect (ITEP-OPE) was addressed via construction of the analytic coupling a{Q'^) = 
at'^'\Q'^)/'K by starting from beta functions /?(a) analytic at a = and performing integration of the corresponding 
renormalization group equation (RGE) in the complex plane. It then turned out that, in order to avoid the 
occurrence of Landau singularities of a(Q^), it was virtually necessary to impose on the coupling a{Q^) analyticity 
at = 0. We tried the construction with many different /3 functions which fulfill such conditions and which, at the 
same time, give relatively tame perturbation renormalization scheme (RSch) coefficients c„ = /3n/Pa (n — 2,3, . . .), 
i.e., where the sequence {|c„|,n = 2,3,...} is not increasing very fast. It turned out that all such beta functions 
resulted either in analytic coupling a{Q'^) which gave r-r < 0.16, significantly below the well-measured experimental 
value r7-(exp.) — 0.203 ± 0.004 of the (strangeless and massless) r,-, or the coupling a{Q^) gave Vr > 0.16 at the price 
of developing Landau singularities. 

This persistent problem was then addressed by a specific modification of the aforementioned beta-functions, 
Eqs. (41)-(42), introducing in /3(a) complex poles and zeros on the imaginary axis of the complex a plane close 
to the origin. In this way, the correct value Vr = 0.203 was reproduced, and the analyticity of a{Q^) and the ITEP- 
OPE condition were maintained. However, the sequence of perturbation RSch coefficients {|c„|, n = 2, 3, . . .} in such 
cases increases very fast starting at n = 4. As a consequence, in such cases the analytic evaluation series of QCD 
observables (including Tt) starts showing strong divergent behavior when terms ~ 05 ~ are included, because the 
coefficients at such terms become large. It remains unclear how to deal properly with this problem. 

In this work we evaluated, in the aforementioned anQCD models, the (timelike) observable r,- and the spacelike 
observable Bjorken polarized sum rule (BjPSR) dBj{Q^) at low Q^, by evaluating only the leading-twist contribution. 



In Ref. |24| an anQCD coupling Ai was constructed directly (not from a function Ansatz) which fulfills the ITEP-OPE condition. 
The construction was performed in a specific RSch and contains several adjustable parameters. Physical observables were not evaluated. 
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and accounting for the chirality- violating higher-twist OPE terms by estimating and subtracting those "mass" terms 
in the case of (see Appendix |b|. This means that the chiraUty-conserving higher- twist contributions, such as the 
gluon condensate contribution, were not taken into account. While the values of the chirality- violating condensates are 
known with relatively high degree of precision and are expected to be the same in perturbative QCD (pQCD-hOPE) 
and in anQCD (anQCD-f OPE), the values of the chirality-conserving condensates have in pQCD-fOPE very high 
levels of uncertainty. For example, the dimension-four gluon condensate, which is the numerically relevant chirality- 
conserving condensate with the lowest dimension in the evaluation of r^, acquires (in pQCD-fOPE) value almost 
compatible with zero: (aG^j,) = 0.005 ±0.004 GeV'^ [57], obtained by fitting pQCD-|-OPE evaluations of the current- 
current polarization operators with the corresponding integrals of the experimentally measured spectral functions of 
the T-decay. In anQCD models, before fitting, the value of (aG^j^) is a free parameter. In principle, the inclusion 
of this parameter, i.e., inclusion of the corresponding dimension-four term in the anQCD-|-OPE evaluation of can 
give us the correc t va l ue o f r^- once the value of the parameter is adjusted accordingly, without the need to perform 
the modification (41 1-(42) of the beta function. It appears that the resulting value of this parameter (aG^^) in such 
anQCD models will be large, especially since it enters the dimension-four term for with an additional suppression 
factor a. Another, more systematic, approach fSS" would be to extract the value of (aG^^), in anQCD models presented 
here, by performing analyses similar to those of Refs. |57[ 158) . involving r-decay spectral functions and suppressing 
the OPE contributions with dimension larger than four by employing specific (finite energy) sum rules. One of the 
attractive features of the anQCD models presented in this work is that most of them give results very similar to each 
other [for a(0), Mthr, ?'t, BjPSR ~ see Tables. iTpV for nonmodified, and V |VII for modified j3 functions] when the 
f{Y) function appearing in the /3 function has various different forms, of the type P[l/0], P[l/1], or EE. 
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Appendix A: Implicit solutions of RGE and singularity structure 



It is evident that for an arbitrary choice of P{F), even when constrained by conditions (21)-(24|, RGE Eq. (15) 
cannot be solved analytically and one has to resort to numerical methods. On the other hand, if one concentrates on 
the question of for which type of /3 function the resulting coupling may have no Landau singularities, more general 
statements can be derived by analytic methods as shown below. 

We suppose that the /3 function has the form Eq. (21) of Sec. |ll] We will show that, if f(Y) of Eq. (21) is any 
rational function (Pade) of type P[M/A^] (with real coefficients and M > TV— 1), with the = analyticity condition 
(24) fulfilled, then there exists in the physical z stripe of F{z) of Fig. fl] (— tt < Imz < tt) at least one pole Zp of F{z) 
\F\Zp) = oo] such that Im(zp) = — tt. The latter means that this is a physically acceptable pole of a{Q^) for < 0, 
i.e., not a Landau pole. The function f(V) being a Pade of the type P[M/N]{Y) means 



f{Y) = f{l/t) = 



(1 - h/t) ■ ■ ■ (1 - tM/t) 

(1 - ui/i) • • • (1 - ujvA) ' 



(Al) 



where the normalization condition /(I) = 1, a consequence of the pQCD condition Eq. (23), is evidently fulfilled 



The fact that this Pade has real coefficients must be reflected in the fact that the zeros tj are either real, or (some of 
them) appear in complex conjugate pairs, the same being valid for the poles Uj. When using the form (Al) in the /3 



function (21 ) and the latter in the integral (27) of the implicit solution of RGE, we end up with the following integral: 



1 



ao/F{z) 



dt t 



M-N+l 



(t - Ml) ■ ■ ■ (t - Mjv) 
{t~h){t-h)---{t~tM) 



(A2) 



where = 1 in the value coming from the first factor (1 — y) in the 13 function Eq. (21). When M > N — 1, the 
integrand in Eq. (A2) can be split into a sum of simple partial fractions l/{t — tj) 



1 



rao/F{z) M 

p / dtll + y^B 

Poao Jao/ai^ I 



1 



(A3) 
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where 



with 



(A4) 



Dn 



- {h-h){h-t2)---{ta'tM) ■ 



(j =G,1,...M) , 
• ■ {tj - tAi) {] 



■M) , 



(A5) 
(A6) 
(A7) 



These formulas can be obtained by direct algebraic manipulations, or by using a symbolic software. Integration in 
Eq. (A3) then gives the following implicit solution of the RGE for F — F{z) in the form z = G{F): 



1 



1 



1 



/3o \F{z) fli, 



M 



3=0 



ao/F{z) - tj 



(A8) 



Within the sum on the rhs of Eq. (A8), the term with j = is (using to = 1) 



ao/flin - 1 



(l-<i)---(l-iM) 



(A9) 



Comparing Bq with f{Y) in Eq. (All we realize that Bq = 1//(1). Consequently, the — analyticity condition 
(24) yields Bq — fSoag [where oq = a{Q^ = 0)]. Therefore, the total coefficient at the j — logarithm on the rhs of 
Eq. ( A8 1 is equal exactly to 1 



1 



Bo = 1 



(AlO) 



On the other hand, this implies that the pole locations Zp at which F{zp) = oo are given by 



z„ = i - 



A ail 



+ ln(-l) - In 



- 1 



l3oao 



(All) 



Let us now investigate where these poles can be localized in the z-plane. In the cases considered here, we have 
< flin < ao [= a{Q^ = 0)], because otherwise (i.e., if < ao < flin) the resulting coupling would give significantly too 
low values of low-energy QCD observables such as the semihadronic r decay ratio"'^^ (tt) or the Bjorken polarized sum 
rule (BjPSR) at low positive Q^'s. Therefore, Oo/oin > 1- In the following, we discuss several scenarios for locations 
of poles Zp-. 



1. If, on the one hand, the roots tj are all real negative, then in the sum over j's (j > 1) on the rhs of Eq. (All ) 



all logarithms ln[— ij/(ao/ai„ — tj)] are un ique and real, as are the coefficients Bj. Hence, this sum is real. The 



only nonreal term on the rhs of Eq. (All ) is ln(— 1) — —in + z27rn. Therefore, Imzp — — vr. This means that in 
such a case there is only one pole and this pole lies on the timelike Q^-axis {Q^ < 0); hence, no Landau poles. 
One of such cases is the one illustrated in Figj2j^a) of Sec.|llj i.e., the case of f(V) being P[l/0] (r2 = 0; M = 1, 
iV = 0) with ti w -1.338. 

2. If, on the other hand, some of the roots tj appear as complex conjugate pairs, the sum over j's (j > 1) on 



the rhs of Eq. (All ) can be real and the same conclusion would apply. However, that sum can turn out to be 



nonreal and we end up with Landau poles. How can this occur? If, for example, tj+i 

imply = B* . However, the corresponding logarithms for j and j + 1 in the sum on the rhs of Eq. (All ) 



t*, thenEqs. (|A4|)-JA7|) 



^® It can be deduced from Appendix [d| Eq. (D13i and Fig. |l3| there, that Fr(t) < 1 and thus the leading-/3o (LB) contribution to Tt is 
r^'^^' < ao. On the other hand, ain = a((3m,c)^) < 0.075. Hence, when < ao < Oin, we have r^^^ < 0.075, significantly too low to 
achieve r^- ft: 0.20. 

^"^ Note: — 7r < Im^ < vr is the physical considered stripe in the complex z-plane. 
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are not necessarily complex conjugate to each other, but can have a modified relation due to nonuniqueness of 
logarithms of complex arguments 



In 



In 



Oo/ O-in — tj + 1 



i27rn,- 



(A12) 



Here, integers nj can be nonzero, but their values must be such that the requirement is fulfilled so that Zp is 
within the physical stripe: — tt < ImZp < tt. Thus, in this case, we can get several poles, some of them with 
—IT < Imzp < TT, i.e., Landau poles. This case is illustrated in the case of f{Y) being cubic polynomial (P[3/0]) 
in Figs. 12 (a), (b), for the case of two different complex values of roots ti: ti = 1 + iO.5 and ti = 1 + iOA. 
Here, the root t2 is then complex conjugate of ti; and ^3 is determined by the pQCD condition ([23]) and 
turns out to be negative. We can see that in the case ti = 1 + iO.5 there are no Landau poles, just a pole 
at Zp = —11.6312 — in. The numerical test with the use of dispersion relation (38) of Sec. |ll] (cf. also Table 
[l| also confirms that a(Q^) = F{z) is analytic in this case. However, in the case = 1 + iOA there are, 
beside the pole at Zp = —10.5023 — ztt. Landau poles at z = —6.32336 ± i2.6005. This can be understood 



P[3/0]; 



1. + 0.5 f 



P[3/0]: r, = 1. + 0.4 j 




FIG. 12: (a) \l3{F{z))\ as a function of z = x + for the beta-function l |21| with f{Y) being cubic polynomial with ti = 1 + iO.5 
[t2 = l- i0.5,t3 = -3.67591); (b) the same as in (a), but with ti = 1 + iOA {h = 1 - iOA,t3 = -3.98969). 



in the following way. The expression for the location of poles Zp is given by Eq. (All I, with the sum there 



over j — 1,2,3. Usually softwares such as MATHEMATICA give for logarithms InU of complex arguments 



U expressions with imaginary part —tt < Im(lnJ7) < tt. In this case, if only the term ln(— 1) in Eq. (All) 
gets replaced by [ln(— 1) — i27r] — —in, the resulting Zp has Imzp = —in, in both cases ii = 1 -I- iO.5 and 
ti = 1 + iOA. Namely, Zp = —11.6312 — in and Zp = —10.5023 — in, respectively. However, if we, in addition, 
replace ln[— ^2/(00/0111 — ^2)] by ln[— f2/(ao/ain — ^2)] + i27r, we get in the case oi ti = 1 + iOA a pole location Zp 
inside the physical stripe —n < Imz < n: Zp = —6.32336 — i2.6005, which is the location of one of the Landau 



poles seen in Fig. 12 ^b); the other Landau pole is at Zp — 6.32336 + z2.6005. 



In general, by adding to each of the logarithms of complex arguments in Eq. ( All[ ) multiples of i2n, we end up 



with a set of possible pole locations Zp. Only those values which lie within the physical stripe —n < Imz < n 
are candidates for the location of (Landau) poles. However, in practice, only some of them represent poles 
F{zp) = 00, while others may have finite values of F{zp). This is so because the RGE integration, for the 
physical stripe of z's, with a specific initial condition aX z = 0, will not cover all the possibilities of these 
multiples. 

Yet another possibility is to have some roots tj real positive. Since we have ao = a{Q^ = 0) by our notation, 
the value a = oq is a root of the beta function 13(a), and there are no other roots of 13(a) in the positive interval 
< a < Qq [note that (3(0) = by asymptotic freedom]. Therefore, we are not allowed to have tj > 1 since 
this would imply that aj — a^/tj < oq is a root of 13(a); hence if tj is positive it must lie in the interval 
Q < tj < 1. Such ij's then fulfill the relations (Q < tj < \ < ao/flin) and hence give a nonreal value of the 



logarithm ln(— tj/(ao/ain — tj)) in Eq. (All I; the value of Bj is real. Therefore, in such a case we generally 
obtain ImZp 7^ — tt, i.e., we generally obtain a Landau pole. 

4. We may obtain Landau poles, or Landau singularities, in several other cases, e.g., when some of the poles u^. 
of the beta function are larger than unity. However, a systematic (semi-)analytic analysis of these problems 
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appears to be too difficult here. We just mention, as an aside, that the appearance of Landau singularities [e.g., 
finite discontinuities of F(z)] usually implies the appearance of Landau poles [infinities of F(z)]. 



When M < N — 2, the implicit solution of the type (|A8|) obtains additional terms on the rhs 



- 3) [if M = - 2: only In F{z] 



In this case the poles |i^(2:p)| 



HF{z)), F{z), 
oo are reached at Zp = — cx). 



0. This implies that in such cases the condition a{Q = 0) = oq < oo cannot be fulfilled. 



Appendix B: Massless part of the strangeless tau decay ratio 

At present, the most precisely measured low-energy observable referring to an inclusive process is the ratio Rt{A.S = 
0), which is proportional to the branching ratio of r-decays into nonstrange hadrons. Consequently, it plays a central 
role for testing the validity of our auQCD approach. However, for a a careful comparison of the available experimental 
result with our theoretical prediction it is essential to extract from the quantity Rr{AS — 0) the pure massless QCD- 
canonic part = rr{AS = 0, uiq = 0). This analysis has already been presented in Appendix E of Ref. fl4'. Here we 
redo it, but with updated experimental values of Rr{AS = 0), of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
element \Vud\ and of higher-twist contributions. The strangeless (V-l-A)-decay ratio extracted from measurements by 
the ALEPH Collaboration [SH [5S] and updated in Ref. [SS] is 

RriAS^O) . ^i:: ^ --h^drons(,)) _ 

= 3.479 ±0.011 . (B2) 

The canonic massless quantity rT-(A5' = 0,TOg = 0) is obtained from the above quantity by removing the non-QCD 
[CKM and electroweak (EW)] factors and contributions, as well as chirality-violating (quark mass) contributions 

rr{AS = 0, m, = 0) = ^ , - (1 + S'^^) ~ Sr^AS = 0, m„,, ^ 0) . (B3) 

o\ Vud\ + (JEWj 

This quantity is massless QCD-canonic, i.e., its pQCD expansion is rr{AS — 0, = 0)pt = a + 0{a^). The updated 
value of the CKM matrix element \Vud\ is [39] 

IKdl = 0.97418 ± 0.00027 . (B4) 

The EW correction parameters are 1 + Sew = 1.0198 ± 0.0006 [51 [55] and S^^ = 0.0010 [53]. The (V+A)-channel 

corrections Srr{AS — 0, m„_d ^ 0) due to the nonzero quark masses are [SHIISS] the sum of corrections (^i^V+'^i^ a)/^ 
with dimensions D = 2,4,6, and 8. It appears that, among the chirality-nonviolating D>2 contributions, the only 
possibly nonnegligible [S^ is the D = A contribution Si^qq-^ = (11/4)q!^ (m^)(aGG)/m^ from gluon condensate. The 
authors of Ref. [55] obtained from their fit the gluon condensate value {aGG) — (—1.5 ± 0.3) x 10^^ GeV*, giving 
thus ^(GG> ~ ~5 X 10"''; their entire value of higher dimension contributions (2 < Z? < 8) to rr{AS = Q,mu^d 7^ 0) 
is (—6.3 ± 1.4) X 10~^. On the other hand, the value of the gluon condensate may be compatible with zero; e.g., the 
T-decay analysis of Ref. [57] based on sum rules gives {aGG) = (0.005 ± 0.004) GeV^ which is almost compatible with 
zero. In our analysis we assume that this is the case, i.e., zero value of the gluon condensate. With this assumption, 
the higher dimension contributions to r-riAS = 0, niu^d 7^ 0) are only the chirality-violating (i.e., due to nonzero quark 
mass) terms, their value being thus 

drr{AS = 0, rriuA 7^ 0) = (-5.8 ± 1.4) x 10"^ . (B5) 



Using the aforementioned results in Eq. ( B3 ) leads to 

r^iAS = 0, = 0)cxp. = 0.203 ± 0.004 , (B6) 
where the experimental uncertainties were added in quadrature. The uncertainty here is dominated by the experimen- 



tal uncertainty SRr — ±0.011, Eq. (B2|. The central value (B6| would increase to 0.204 if th e glu on condensate value 
(aGG) = (-1.5 ± 0.3) x 10"^ GeV^ Ref. ^563 was taken. The central value 0.203 of Eq. is also obtained by 



using the analysis and results of Ref. [57 , but with the updated values Rt-{AS — 0) of Eq. (B2) and \Vud\ of Eq. (B4| 
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Appendix C: Higher order terms in analytic QCD 

Here we summarize the general approach to calculate higher order corrections in analytic QCD (anQCD) models, 
as described first in our earlier works [131 [14]. In order not to confuse the general analytic coupling a{Q^) with pQCD 
coupling apt(Q^), we will use in this Appendix the notation Ai{Q^) for the analytic coupling. 

First we note that the analytic coupling Ai{Q^) does not fulfill the ITEP-OPE conditions ^ in any of the anQCD 
models that have appeared in the literature up to now.^® Nonfulfillment of ITEP-OPE conditions implies that the 
respective beta function l3{Ai) = dAi{Q^)/dlaQ^ is not analytic in Ai (cf. arguments in Sec. [n]). Consequently, in 
these models the beta function, which is usually not known explicitly, cannot be Taylor expanded around Ai = 0, and 
therefore the powers Ai cannot be expected to be the analytized analogs of aj^j. In fact, they usually are not. The 
construction of AniQ^), the analytic analogs of apt(Q^)" (n > 2), is yet another important ingredient in anQCD. 

A spacelike massless observable V{Q^), in its canonical form, has the following perturbation series: 

1^{Q\t = apt + dial, + ds^t + ■ ■ • , (CI) 
and the corresponding truncated perturbation series (TPS) is 

2?(Q2)[^1 = apt + dial, + ■■■ d^-ia^, . (C2) 

Here, Opt and d^'s have given renormalization scale (RScl) and scheme (RSch) dependences. Analytization means, 
in the first instance, to replace in the first term apt by Ai{Q^). For treating the higher order terms, there are, in 
principle, several options at hand. For instance, one could replace all powers of Opt by the corresponding powers 
of Ai {apt H' Ai). Or, as is done in MA, one could subject each a^, to an analogous analytization procedure as 
Ai (if such an analogous procedure unambiguously exists), yielding additional analytic couplings Qp, i— >■ An, where, 
in general. An 7^ Ai- In MA such a prescription unambiguously exists. The advantage of such a prescription in 

MA lies in the fact that the RGEs governing the running of a1,"^)'s, as well as the RSch dependence of ^n'^^^'s, 
are identical to the corresponding pQCD RGEs and RSch dependence once the replacements ap, 1— >■ -AI'^^^ are 
performed there j64| . We consider this property as physically important, especially because there is a clear hierarchy 
^(MA) ^ |^(MA)| ^ |^(MA)| _ positive values. Among other things, this hierarchy implies that the 



MA-analytized version of the TPS Eq. (C2) 



m'iZ) - -^r^^ + ^i4''^^ + • • • d^-iAT^ , (C3) 
becomes systematically more RScl and RSch independent when the truncation index N increases 

ai?(Q2;RS)[^l 



d{RS) 



= fc^-4^Ti +0«r2)- (C4) 



Here, "RS" stands for logarithm In/i^ of RScl /i, or for any RSch parameter Cj = fij/Po {j > 2). 

However, when constructing anQCD models beyond MA, by changing the discontinuity function piicr) = 

Imapt(— cr — ie) appearing in the dispersion relation ^ for Af^'^\Q'^) [TTJ [T31 [H], or by different constructions 
of Ai{Q'^) (cf. [7l-[T0l [T2] and references therein), the meaning of "analogous analytization" of higher powers a^t be- 



comes unclear or, at best, ambiguous. On the other hand, it is almost imperative to maintain relations (C4) in any 
anQCD model with hierarchy Ai > \A2\ > l^sl • • • , because then the physical condition of RScl and RSch indepen- 
dence of the evaluated observables is guaranteed to be increasingly well fulfilled at any when the number of terms 
increases. 

Furthermore, it is preferable to have the higher power analogs a^, An not simply constructed as An = (^1)", 
but rather by application of linear (in ^1) operations on Ai, such as, e.g., derivatives and linear combinations thereof. 
The underlying reason is the compatibility with linear integral transformations (such as Fourier and Laplace) [65) . In 
linear transformations, the image of a power of a function is not the power of the image of the function. 



Except for Ref. |31 | where some of the main results of the present work have already been summarized, and Ref. |24| where a direct 
construction of an analytic coupling Ai with several parameters was performed (cf. footnote |15| in this work). The anQCD model of 
Ref. |12l fulfills this condition approximately. 

Such a construction of AniQ^), as a linear operation applied on J[i{Q^), was presented in anQCD in Refs. I13II14||20) . 
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The construction of higher order analogs An (appUcable to any anQCD model) which obey all these conditions was 
first presented in Refs. [13l[T4]. The procedure proposed there for obtaining An from a given anQCD coupling Ai, in 
a given RSch, is the following: First we define the logarithmic derivatives of Aii^i^) (where — kQ^ is any chosen 
RScl), i.e., we define 



An+l{fl^) (= a„+i(/x^)) 



(-1)" 

I3^n\ 9(ln^2) 



2\n 



(n = l,2,...) 



(C5) 



In order to understand the following construction of AnS given below, it is convenient to consider first the corre- 
sponding logarithmic derivatives in pQCD 



( 2^ (-l)"a"apt(M^) 



(n = l,2,...) 



These^^ are related to the powers 0"^ via relations involving the Cj coefficients of the pQCD RGE Eq. (14 



Opt, 2 = Opt + ciflpt + C2apt H , 

Opt, 3 = apt + ^ ' 

apt.4 = Opt H , etc. 

The above relations are obtained by (repeatedly) applying the pQCD RGE. The inverse relations are 

.'5 

'^pt — ^pt,2 — Ciapt,3 



"1 - C2 apt,4 



'^pt — ^pt,3 — ^Cittpt^i 
^pt — Q^pt,4 + • • • , 



etc. 



Now we adopt the following replacement on the rhs of Eqs. (Cf0l-(C12): 

Opt Ai , apt,„+i An+i (n = f , 2, . . .) , 
and use the generated expressions as definitions of An, the higher power analogs of pQCD powers a^^ 

'5 



A2 ^ A2- C1A3 + ( - C2 ) A4 



■^3 = -^3 - -C1A4 

Ai ^ Ai^ , 



etc. 



(C6) 

(C7) 
(C8) 
(C9) 

(CIO) 

(Cll) 
(C12) 

(C13) 

(C14) 

(C15) 
(C16) 



It is then straightforward to see that the analytic ("an") series obtained from the perturbation series (CI) via 
replacements apt 1— Ai, a^^ 1— )■ An 



2?(Q^)an = ^1 + diA2 + d^A^ + • • • , 
gives the corresponding truncated analytic series 

I?(Q2)[^] =Ai + diA2 + ■■■ dN-iAN , 



which really fulfills the condition (|C4|) of increasingly good RS-independence, now in any anQCD model 

kpfAN+i + 0{An+2) , (RS = In^^;c2;c3; . . .) . 



9P(g2;RS)Ln' 



5(RS) 



(C17) 
(C18) 

(C19) 



An expansion of the Adler function in terms of apt,„+i(/i^) is used in Ref. 66 for an evaluation of Tt in the context of pQCD; this 
"modified" contour improved perturbation theory (mCIPT) was shown there to have advantages over the standard (CIPT) approach, 
most notably a lower RScl dependence of the result. 
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This relation continues to hold even if we truncate relations (C14)-(C16| at the order ~ An (including the latter). 

The above presentation suggests that, instead of the perturbation series (Cll in powers of Opt, a modified pertur- 
bation series in logarithmic derivatives apt,n+i (C6) can be used 



whose truncated form is 



2'(Q^)mpt = Opt + diapt,2 + d2apt,3 



V{Q') 



nipt 



where "m" in the subscript stands for "modified," and the modified coefficients dj [j = 1, 
the original coefficients dj 



di 
d2 

ds 



di 
d2 

ds 



cidi , 



C2 Ml , 



etc. 



(C20) 
(C21) 

, — 1) are related to 

(C22) 
(C23) 

(C24) 



When applying analytization to the modified perturbation series (C20), via replacements (Cf3), we obtain modified 
analytic series ("man") 



whose truncated version is 
Its RS dependence is 



[N] 
man 



2?(Q^)ma„ = ^1 + diA2 + d2A3 + 



T^iQ^t'Hn =-^1+ dlA2 + ■■■ dN-.AN 



kNAN+i+0{AN+2) {-^ An+i) , (RS = In^^;c2;c3;...) 



(C25) 
(C26) 

(C27) 



a(RS) 

It is interesting that in virtually all anQCD models [i. e., mo dels that define yli((5^)] holds the hierarchy Ai > \A2\ > 
lAsl > ■ ■ ■ at (almost) all complex Q^. Therefore, Eq. (C27) signals an increasingly weak RS dependence of 'D{Q^)[^}ri 
when N increases, at any value of and RScl 

We stress that the analytic ("an") and modified analytic ("man") series [Eqs. (C17) and (C25l, respectively], if 
they converge, are identical to each other due to relations (C22)-(C24) and ( |C14 )-(C16). 

In the specific case of M A, i.e ., whe n Ai = .4^^^'' of Ref. [4J, it can be shown (using the results of Ref. [64j) that 

the above procedure, Eqs. (C14)-(C16), gives the same higher power analogs A^^^ as the analytization procedure of 
Ref. '5 (APT) that uses the MA- type dispersion relation involving Imapj((3^ = —a — ie) 



1 /-oo (pt) 
Jo + 



(^) 



(C28) 



where p'n^\a) — Imapt(— ct — ie) (n = 2, 



We note that -41^'^'' ^ {Af''^')"'- Furthermore, construction of An 

„(pt) 



,(MA) 



according to relations ( C14 1-( C16 ) in other models of anQCD (e.g., where Ai is constructed from a modified pi =^ 
e.g. Refs. [TI] [131 [H] ) also in general leads to An ^ -4". However, if analytic Ai{Q^) = a{Q^) is constructed from 
RGE with beta function /3(a) analytic at a = 0, as is the case in the present work and Ref. ^31], it is straightforward 
to see that construction ( C14[ )-(C16) gives 



Ar, 



(n = l,2,...) . 



(C29) 



In those anQCD models of analytic Ai{Q'^) where the aforedescribed construction gives An 7^ A'l for n > 2 (such 
models do not appear in the present work), using yl" instead of An is not a good idea for at least two reasons: (1) 
such a construction is formally not linear in Ai [see the discussion before Eq. (C5|]; (2) the RS dependence of the 
resulting truncated "power" analytic series 



■^pan (Q 



2^[N] 



Ai + diAl 



■ d]s-iAi 



N 



(C30) 
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is not entirely analogous to Eq. ( |C19 ) or Eq. (C27), but is rather 



NP 



(AT) 



(C31) 



where NP(jv) is an increasingly complicated expression of nonperturbative terms (such as when N increases, 

and |NP(7v)| in general does not decrease when N increases. 

Appendix D: Leading-/3o (skeleton-motivated) resummation in anQCD 

First we summarize here the resummation formalism for the leading-/3o (LB) part of inclusive spacelike QCD 
observables in anQCD models, as presented in [13 [H]. Subsequently, we present application of this formalism to LB 
resummation for the Bjorken polarized sum rule (BjPSR) dBj(Q^) and, in a newly modified form, to the r decay ratio 



Massless spacelike QCD observables 'D{Q-), in canonical form, have the pQCD ("pt") expansion (CI I in powers of 
Opt, where apt = apt(/i^; C2, . . .) is defined at a given RScl ^ and in a given RSch (c2, C3, . . . ). In the scaling definition 
of jjL we use the convention A = A, which is the MS reference scale for RScl's /i [the so-called V scheme Ay is related 

2 

to A via A = Ayexp(C), where C = —5/3]. The considered RSch classes will be such that the RSch coefficients 
Pk = PoCk {k > 2) are polynomials in n/, and consequently in /3o — (11 — 2n//3)/4 



Pk = PoCk 



bk.Pl (A: ==2,3,...) 



(Dl) 



We recall that /3o — (11 — 2n//3)/4 and /3i = (102 — 38n//3)/16 are both universal (RSch-independent) parameters. 
RSch's MS and 't Hooft are clearly special cases of such RSch's. The RSch i ndep endence of 'D{Q^) implies a specific 
dependence of coefficients c?„ on the RSch parameters jl?]; this and relations (Dl ) imply that the coefficients (i„ have 
specific expansions in powers of /3o 



.(1) 



(D2) 



k=~l 



We note that c'l'Li = 0. In MS RSch, the negative power term cx l//9o does not appear. Relations (D2) and (C22)- 
(C24l) imply that the modified perturbation ("mp t") e xpansion (C20l of 'D{Q^) in logarithmic derivatives Opt.n+i of 



( |C24|) im 
Eq. feel 



have coefficients d„ of a form similar to ( D2 ) 



E 

fe=-i 



^nk PO 



Specifically, the leading-/3o terms in Eqs. (D2) and (D3) coincide^^ 



^1) 



.(1) 



(D3) 



(D4) 



The LB resummation of the inclusive spacelike I?((5^) is obtained in pQCD via integration of apt(/i^) over various 
scales — tQ^ exp{C) and weighted with a characteristic function Fj,{t) according^^ to formalism of Ref. ^67 



(lb; 
pt 



Jo 



dt 



(D5) 



The integration cannot be performed unambiguously, due to the Landau poles of apt at low values of t. In anQCD 
Opt here is simply replaced by analytic Ai (= a) 



dt 



(D6) 



21 Note that /3i = bio + 6ii/3o (with: feio = -107/16 and 611 = 19/4); therefore, ci = /Bi/I^o is ~ /3g in the leading-/3o (LB) hmit. 
The superscript S indicates here that the observable is Euclidean, i.e., spacelike. 
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where now the integration is unambiguous since there are no Landau poles. Expansion of the analytic coupling 
Ai{tQ'^e^) around the RScl scale /i^, i.e., Taylor expansion in powers of £ = ln[/x^/(iQ^e'')], gives 



(D7) 



n=l 



We thus see that integral (D6), in anQCD, represents exactly the leading-/3o (LB) part of the modified analytic 
( "man" ) expansion ( C25 ) in Appendix [C] The truncated series of the latter is given in Eq. ( C26 1 . We stress that 
the above expansion is performed at a given RScl /i and in a given RSch [02,03, . . . - cf. Eq. (Dl)]. In anQCD it is 



convenient to perform explicitly the LB resummation (D6) since the integral there is finite, unambiguous, and RScl 
independent. 

The characteristic function F^{t) for the BjPSR 'D{Q'^) — d'Q^{Q^) was calculated and used in Ref. [13] (on the 



basis of the known ;68j coefficients cl^2 for it), and was presented in Ref. [13] 



(8/9)t [1 - (5/8)t] t<l 
(4/(9t)) [1 - l/(4t)] t>l 



(D8) 



The (nonstrange massless) canonical^"^ semihadronic r decay ratio r^- = r7.(AS' — 0,mq = 0) is a timelike quantity, 
and can be expressed in terms of the massless current-current correlation function (V-V or A-A, both equal since 
massless) [55] 



ds 



imn(g2 



(D9) 



Use of the Cauchy theorem in the plane and then integration by parts leads to the following contour integral form 

MM- 



1 

2tt 



(l + e^0)3(i_e#) dA^iiQ^ = my*') 



(DIO) 



with dAdii Q^) = —dIl{Q^) /d\nQ^ being the massless Adler function. In pQCD, use of the Cauchy theorem to the 



expression (D9| is formally not allowed. This is so because npt(Q^), being a power series in apt(Q ^) [or: apt(K(5^)], 



has Landau singularities along the positive axis < < A^. In pQCD, expressions (D9) and (DIO I are two different 
quantities; in anQCD models they are always the same. 

The massless Adler function (iAdi(Q^) is a spacelike (quasi)observable. On the basis of the known coefficients Cnn 



for it [70] [7T], its characteristic function FAdi(T) was obtained in Ref. |57], and from it and using relation (DlOl the 
characteristic function for was obtained in Ref. [72], in the timelike LB form 



-(AS" = 0,m^ = 0) 



(LB) 



dt 



(Dll) 



Here, the superscript A4 indicates that these are Minkowskian (timelike) quantities; 2li is the timelike coupling 



2li(s) 



da 



Pi(ct) ; 



(D12) 



and the characteristic function F^{t) was obtained in [72].^^ 

It turns out that, in the calculations in the present work, it is inconvenient to calculate the LB-contribution to 
by using formula (Dll) which involves function 2li(s). This inconvenience consists in the following: in this work, 
RGE ( [l5| Eqs. (16)-(17l] is integrated in the entire physical stripe in the complex z plane, and as a result of this 
we numerically obtain, among other things, the quantity pi(cr) = Ima((5^ = —a — ie) — lmF[z = \z\ — in) [with: 
\z\ = ln{a / fifjj]; to obtain the quantity 2ti(s), yet another numerical integration (D12) is needed, and then we go 



Canonical form, in the sense that its pQCD expansion is = apt + C'(apj). 

In fact, the quantity Wt of Ref. |72| is related to here via: F^{t) = (t/4)WT(i). Full expression for F^{t) is given in Eqs. (ClO)- 
(Cll) of Ref. |14| : however, a typo appears in the last line of Eq. (Cll) there: in a parenthesis there, the term +3 should be written as 
3t^; the correct expression was used in calculations in Refs. |13l I14| . 
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with this 2li(s) into the integration (Dll 



precision of calculation is expected to be 



There are too many successive numerical integrations involved, and the 



Therefore, we perform in integral (Dll) integration by parts, using relation (M.i{s)/d\iis = —pi{s)/7r [cf. Eq. (D12)], 

(LSI 



and we obtain the expression of ri- in terms of the discontinuity function pi{s) 

1 f°° dt 



,(LB) 



t 



Frit) piite'-mi) , 



where 



Frit) 



* dt' 



I -vrFnt') 

^ 



(D13) 



(D14) 



Integration in (D14| can be performed analytically, and the result for Frit) is iCp — 4/3): 



Frit)/iACF) = -^Li2(-<) (t'^ + et^ + ISt^ + 10<- 12nn(<) - 3) - 2tLi3(-i) 



17^ { - 72 ln(<) [t (-2<2 - 47t + 6) + 2 (t^ + 6<^ + ISt^ + lOt - 3) ln(< + 1)] 
259i'' - eOOt^ - 6948^2 - 5184iC(3) + 7344t + 72(t + &)t^ In^(t)} (t < 1), 



(D15) 



Frit)/i^Cf 



1 

432 



mit^ + Qt^ + m - 2)Li2 ( -7 ) - 108Li2(-t) + 864tLi3 



+432tLi2 - 1) - 9 [(t^ + 8t + 36)i2 + 96] ln(i + I) - t [9t(4t + 23) + 598] 

-18 [2tit^ + Qt^ + 18< + 22) - 3] \T?it) 

+3 [(3t^ + 12^3 + A2t^ - IMt + 111) + 12(t2 + 4< + 9)(t + 1)^ ln(i + 1)] ln(t) 

1 



+9t(t^ + 8t^ + m - 96) In - + 1 



t 



+432(ln(t) - 2) [t ln(t) - (i + 1) ln(i + 1)] + 648C(3) - IMtt^ + 841 
3C(3) , 2.^ 463 ^ 



(D16) 
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The function Frit) is continuous and monotonously increases when t increases. Its value is zero a,t t — Q, and one at 



t = +00. It is depicted in Figs. 13 as a function of t and Int. 
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FIG. 13: Characteristic function Fr{t) which appears in the LB integral (D13l of r,-: (a) as a function of t; (b) as a function of Int. 
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Appendix E: Inclusion of beyond-the-leading-/3o (bLB) terms in anQCD 



In pQCD, perturbation expansion of any m assle ss spacelike observable T>{Q^) can be written in the form (CI) or 
(C20). In the considered (large) RSch classes (Dl), the coefficients d„ and c?„ can be written in the form (D2| and 



(D3|, respectively. Leading- /^q (LB) resummation (D6| reproduces one part of these terms, Eq. (D7|. In practice, for 
inclu sive spacelike observables only the leading- /3o parts CnnPo of coefficients dn and d„ are known for all n [cf. also 
Eq. (D4|], while the coefficients known in their entirety are only the first two or three: di,d2,d^ di,d2,d^^ 



cf. Eqs. ( |C22l )-( [C24l )]. For this reason, t he m ost that one can include in the evaluation of any such observable in 
anQCD are all the LB contributions, Eq. (D7), and the beyond-the-leading-/3o (bLB) terms of order a^,a^, anda'* (<^ 
of order 02,03,04]. 

In practice, the coefficients ^1,(^2,^3 and 
[c2(MS), C3(MS), . . .] and with^^ RScl fj? = Q^; we will denote such quantities with the bar over them. In gen- 
eral, the evaluations are performed in another RSch (c2, C3, . . .) (e.g., in the present work the RSch as dictated by the 
chosen 13 function used), and another RScl 



Cnn^Q are calculated and given in the literature in the MS RSch 



- exp(C) (C - 1) 



(El) 



The LB contribution (D6) is RScl independent; however, it depends on the RSch. The truncated bLB contribution 



still has some remnant RScl dependence due to truncation, and is RSch dependent. 

The dependence of the coefficients dj on RScl and RSch can be deduced systematically, by the requirement of RScl 
and RSch independence of the observable V and using the known RScl and RSch dependence of the pQCD coupling 
apt(C; C2, C3, . . .) [H]. The resulting dependence of dj is 



dl = dl +/3oC (=di) , 

2/3oCZ+PoC^] -(C2-C2) 



d2 = d2 
ds = ds 



3/3oCd2 + 3(3'^C'''di + /S^C^ 



-3(rfi + /3oC) + -ci 



(C2-C2)- ^(cs-Cg) 



(E2) 
(E3) 

(E4) 



etc. On the other hand, the RScl independence of LB contribution (D6|-(D7) implies for the LB coefficients (D4| the 
following RScl dependence (they are RSch independent) 



k=0 ^ 



fcfc' 



(E5) 



where Cqq' = 1 by definition. When we subtract from the modified analytic ( "man" ) series ( C25 1 the LB contribution 



(D7|, we obtain the bLB contribution separately 



rif — ~ 



t 



(E6) 



where C — —5/3 as mentioned earlier in Appendix [pj An+i are in RSch (c2, C3, . . .) and at RScl — exp(C), and 
the coefficients {Txi)n are 



(E7) 



whe re dn and C nl a re related with the corresponding (b ar) qua.ntitie s in MS RSch and RScl /i^ = via relations 
(p2|-([E4l) and (IesI. This, and apphcation of relations (|C22l)-(|C24| in MS RSch and RScl = Q^, allows us to 



obtain the first three coefficients {Tx))n by knowing the first three coefficients d„ (n — 1, 2, 3) (all c"^^ are known) 



Sometimes, c]ln^ a-re calculated and given in the literature at RScl fi? = exp(C) = exp(— 5/3). 
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Another variant of evaluation of D in anQCD is not to perform the LB resummation (D6) in (E6), but rather use 



its expanded form (D7|. This leads to 



n+1 , 



(E8) 



where a = Ai= Ai{Q'^ ex p(C) ; C2, . . .). Series (E8) was obtained in Appendix [c] in Eq. (C25|. 

In principle, both series (|E6[) and (E8 1 must lead to the same result if the series are convergent. However, in practice, 



only the first three terms in the sums there {n — 1,2,3) are known. Hence the series (E6| and (E8) truncated at 
n — 3 



[4] 



dt 



F%{t) Ai{tQ^e^) + Y,{Tv)rAn+i 



n=l 



2?ma„((3')'^^ = Al+J2dnAn+l , 



(E9) 
(ElO) 



n=l 



will give in general somewhat different results, the difference being ^ .45 As). In theory, the LB-resummed 
truncated version (E9) is better since it includes more contributions than the simple truncated version (ElO). Which 
of the two is better in practice, in the case of a specific considered inclusive observable 'D{Q^), can be decided 
numerically, e.g., by establishing which of the two truncated series has weaker variation under the variation of the 
RScl (<S4> under the variation of C). If is not an inclusive observable ( e.g., jet observables, etc.), LB resummation 

cannot be performed since F^{ t) d oes not exist, and only the expression (ElO) is applicable in such a case. 

The bLB part o f exp re ssion (E6), and the sum over An+i in Eq. (E8), can be reorganized into s ums over .4„+i's 
as defined in Eqs. (C14)-(C16) pT^+i — a"+^ in our paper since /?(a) is analytic in a = 0, Eq. (C29)]. In such a case, 
the truncated analytic expressions analogous to (E9)-(E10) are 



2?ii;B+''LB)(g2)[4] 



dt 



F^{t) Ai{tQ 



3 

E 

n=l 



{iv)nA^ 



■n+1 



(Ell) 



2?an(Q 



2^[4] 



Ai 



(E12) 



The t runcated series (E12) was obtained in Appendix [C| i n E q. ( C18 ). Aga in , theo retically, the truncated expansion 
|e11|) is better than (|E12|). All the truncated expansions (|e9I, (IEIOI), (|e11|, |e12| differ from each other by ~ A 



(1) 



.(i)i 



and Eq. (E7) 



.45. Our numerically preferred version of evaluation will be the truncated expansion (E9). 

Expressions for bLB coefficients (Tx))„ {n = 1,2,3), appea ring i n Eqs. ( |E7[ ) and (E9 ), a re o b tain ed from the (us ually 

known) coefficients (j = 1,2,3) via successive use of Eqs. (C22)-(C24) \dj i-> d^]; Eqs. (E2)-(E4) [d^ ^ d^]; Eq. (E5) 



It turns out that these coefficients are equal to the coefficients i„+i as derived i n Ap pendix A of Ref. [H], tn+i — 
(Tx,)„, as it should be.^^ The bLB coefficients (tx))« ("- = 1, 2, 3) appearing in Eq. ( |Ell| , on the other hand, turn out 

to be equal to expressions t„+i = t^l^^ + ■ ■ ■ ^^"i^l^'' of Appendix A of Ref. [14] when the RScl parameters Ck there are 
all set equal to C. 

In our evaluations of BjPSR and r^., we will use dn [n ~ 1,2,3) coefficients (in MS RSch with RScl /i^ — Q^) for 
massless BjPSR 'D{Q^) — dBj(Q^) and massless Adler function 'D{Q'^) — d^diiQ^)- 



In Eq. (A18) for 44 = {Td)^ of Ref. |14j there is a typo: in the first line the last term should be — <5fe2i3(cJ^' +C) instead of — 562i3Cj'^'. 
The correct formula was used in the calculations there; Eqs. {89)-(92) in Ref. 1141 . which follow from Eq. (AIS) there, are correct. In 
terms of the quantities of Ref. |14| . Eq. (A18) there (without the typo) can be rewritten in the form: 

U = {T-Dh =U- (1/2)(C3 - C3) - (C2 - C2) [scJV + 3(5^1' + C)/3o - (5/2)ci] + 3C/3o?3 + 3C'l3ic[o ■ (^13) 
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Coefficients di and d2 for massless BjPSR were obtained in Ref. |73) . 

(dBj)i = -^+2/3o , 



iJ2 



-35.7644 + 10.5048/3o + 6.38889^^ 



(E14) 
(E15) 



and ds was estimated in Ref. [74] 



(rfBj)3«130 (n/=3). (E16) 

The leading-/3o coefficients c™ for BjPSR were calculated in Ref. [68j in the MS RSch and at RScl = exp(C) 
(where: C = —5/3). When changing RScl to /i^ = using an "inve rted" version of relations (E5| (with cjn i— » ci\|, 



cifc ^ 4fc ' and C -C = +5/3), we obtain c\\' = 2 [cf. Eq. (lE14f]; c^j^ = 115/18(w 6.38889) [cf. Eq. (|E15[)]; and 



7.(1) 



^33 



605/27(w 22.4074). 



Coefficients dn {n — 1, 2, 3) for the massless Adler function were obtained in Refs. [75^77) . respectively 

(dAdi)i = ^ + 0.691772/3o , 

(dAdi)2 = -27.849 + 8. 22612/3o + 3.10345/3^ , 

(dAdi)3 = 32.727 - 115. 199/?o + 49.5237;9^ + 2.18004/3;^ . 



(E17) 

(E18) 
(E19) 



The light-by-light contributions are not included in these coefficients; however, they are zero when = 3, and the 
value Uf = 3 is used in the evaluation of (iAdi(Q^) and subsequently in the evaluation of r,-. The latter observable (with 
AS = and the mass effe cts su btracted) is calculated by using the massless Adler function dAdi(Q^ = m^_exp{i(j))) 
in the contour integration (DlOl. Specifically, applying this contour integration to the analytic expansion (E6| of the 
Adler function, we obtain 



(r.)i^f ™^ = + E (rAdl)„I(^„+l,C) 



(E20) 



where 



/(X+i,C) = 



2tt 



and Vt is given in Eq. (D13). In practical evaluation, the sum in (E20) is truncated at n = 3 

dt 



irr) 



(LB+nLB),[4] 
man 







t 



Fr{t) piite'^ml) + E (rAdi)„/(A.+i,C) 



n=l 



(E21) 



(E22) 



The ot her t hr ee an alytic versi ons of evaluation are obtained by contour-integrating, via (DIO), the analytic truncated 
series (|E10[), (|Ell[) and (|E12[) of massless Adler function V{Q^) = dAd\{Q^)- 



(r V"! 

/ N(LB+nLB),[4] 



/(^i,C) + E (dAdi)„/(A.+i,C) 



dt 



- Fr{t) pi{te^ml) + E {tAd\)nI{An+l,C) 



{Vrtl = /(^1,C)-HE (^Adl)„/(A.+ l,C) 



(E23) 
(E24) 
(E25) 



Again, all four versions of the anQCD evaluation of r^. differ from each other by ~ ~ A^. The truncated expansion 
(E22) is our numerically preferred version. 
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